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SECTION 1 
INTRODUCTION 


The objective of this study was to investigate and develop 
techniques to characterize the quality of the products produced 'by an image 
data system in terms of measurable performance of the various system 
elements. The product quality measures were required to be meaningful 
to the user and related to unambiguous, measurable element performance 
descriptors. A specific objective was to investigate the requirements for 
quality control procedures at the Earth Resources Technology Satellite NASA 
Data Processing Facility (ERTS NDPF) to demonstrate the developed 
approach. 


The fundamental problem was to establish an objective, 
functional measure of image quality and a means to predict that measure 
from accepted system element performance descriptors. Supporting 
tasks undertaken included familiarization with the NDPF, survey of 
ERTS user's "image quality" requirements and survey of currently 
accepted image quality evaluation techniques. It was found that adequate 
performance descriptors existed for the NDPF elements but that no 
adequate measure of image quality existed in terms of user's requirements. 

The measure of image quality developed is the error intro- 
duced by the image data system into estimates of target characteristics 
from measurements made on output products. If the system is characterized 
by a sequence of operators corresponding to element performance descriptors, 
then the parameter estimation errors are mathematically related to the net 
system operator. A software program, Image Data System Simulation 
(IDSS), was written to implement the approach. IDSS uses computer simula- 
tion to synthesize the image data system from its element performance 
descriptors and computes the estimation errors for specific target parameters 
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The IDSS program, applied to the ERTS NDPF, can be used 
to analyze a large class of image processing systems, and should be useful 
in system design or upgrading. 

In addition to NDPF quality control, the ability to monitor 
payload sensor performance by measurements made on NDPF output 
products was examined. Edge gradient spectral analyses software was 
written for NASA's computer as a tool for Modulation Transfer Function 
(MTF) measurement. The feasibility of a technique for correction of 
radiometric distortion due to atmospheric scattering and' sensor effects 
based exclusively on measurements made on output products was also 
considered. 


The next section of this report presents a summary of the 
background information developed at the onset of the study. It includes the 
results of survey of the data flow within the NDPF, ERTS imagery users' 
requirements, and state-of-the-art image quality evaluation techniques. 
Section 3 contains the mathematical bases of the approach and Section 4 
describes the software written for its implementation. Section 5 identifies 
the elements of the NDPF a-nd their associated performance descriptors. 
The results of the application of the developed technique to the NDPF are 
presented in Section 6. Finally, conclusions and recommendations are 
presented in Section 7. 
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SECTION 2 


BACKGROUND INFORMATION 

At the onset of this study a specific operational description 
of the NDPF and its subsystems, particularly the Bulk Processing 
System, had to be developed. In addition, it was necessary to form a 
concept of what the users of image products expect to obtain from them. 
Finally, a brief survey of the state-of-the-art to image evaluation was 
conducted. The results of these surveys provided necessary background 
information that is summarized in this section and presented in more 
detail in Appendix A. 

2. 1 NDPF Description 

This study was limited to consideration of the NDPF video 
data to film product conversion process and digital products excluded. 

The inputs to the NDPF are video tape recordings from the ERTS payload. 
Although payload and telemetry link elements are external to the NDPF, 
the quality of the video tape input will influence the image product quality 
and consequently these elements were represented in the analyses. 

A block diagram of the data flow within the Bulk Image 
Processing system is given in Figures 2-1 and 2-2. This block diagram 
served as a basis for system modeling. The following comments refer 
to the circled numbers in Figures 2-1 and 2-2. More specific data is 
presented in Appendix A. 

The Return Beam Vidicon (RBV) (T) consists of three 
boresighted RCA vidicon systems in three different spectral bands. 
Radiometric and geometric calibration capability exists in the payload. 
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Figure 2-1 NDPF DATA FLOW BLOCK pi AG RAM (PART.1) 
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Figure 2-2 NDPF DATA FLOW BLOCK DIAGRAM (PART 2) 






















The RBV telemetry (2) is analog. The net frequency response 
of the RBV videotape record is given in Figure 2-3 (taken from Ref. 2, 
p. A -5).. 

The Hughes Multispectral Scanner (MSS) ( 5 ) has four con- 
jugate linear detector arrays, each in one spectral band. The MSS tele- 
metry (5) is digital. The frequency response of the MSS videotape record 
is given in Figure 2-4 (taken from Ref. 2, p. A- 13). MSS data is D/A 
converted (*5) in the playback operation. Radiometric corrections are 
applied in the digital domain before conversion. 

The Electron Beam Recorder (EBR) (s) produces all 
archival latent images. The film type used is Kodak SO-43 8 (Ref. 1). 

A fifteen step gray scale is put on each image. 

The first processor ( 9 ) is a Kodak Versamat used only 
for processing archival images. 


The quality control blocks (QC) (10) consist of standard 
Kodak chemistry quality control plus the placing of a special target on the 
head and tail of each roll processed. The target consists of two frames: 
one containing a gray scale, another containing five equal, uniform 
density patches and a standard Air Force tri-bar target. The gray scale 
is read and a Hurte r-Driffield (H-D) curve fit to the data points. The 
density values at two exposure levels are plotted and deviation from 
nominal values used as a processing quality criterion. The constant 
density is read at the five format positions to provide uniformity data. 

The Air Force tri-bar target allows determination of on-axis resolution. 
Additional information on the quality control procedures is included in 
Appendix A: 


6 



©IGINAL PAGE IS {%> a$N0dS3y havm ayvnos 

)F POOR QUALITY] 



4.6 9.2 13.8 18.4 


23.0 27.6 32.2 


36.8 41.4 46.0 (AT 70 mm FILM FORMAT) 


CYCLES/mm 


Figure 2-3 MEASURED RBV SYSTEM SQUARE WAVE RESPONSE 




Figure 2-4 MEASURED MSS SYSTEM SQUARE WAVE RESPONSE 


2.2 Users' Requirements 


It is required that the image quality characterization 
resulting fzom this study contain sufficient. information for users to 
determine the adequacy of ERTS imagery for individual needs. Users' 
ability to define "image quality" is not established and completeness 

of an image quality characterization based on users' requirements would 
not be expected. But a general -understanding, at least, of what tasks 
users of ERTS imagery would like to accomplish is certainly required 
if the "quality"' which is controlled is to have relevance. 

A literature search was consequently conducted. Based 
on the literature sampled, no definition of "image quality" useful in 
accomplishing the study objectives can be drawn from the users. Surely, 
one might adopt the terms "radiometric fidelity", "geometric fidelity", 
"resolution", but these terms have different meanings to different people 
and are certainly not sufficiently well defined to provide a useful basis 
for quality control criteria. They are general terms which classify 
rather than specify the ability to make, certain measurements on 
photographs. 


We examined the users' tasks to determine what sort 
or measurements each user user was making. The results seem to 
span the following questions: 

1. How accurately can a boundary between different 
transmission levels be located on a photograph? 

2. How well can the radiance, size, and location of 
small objects be measured? 

3. How well can the distance between two objects 
or boundaries on a photograph be measured? 

OF POOR OpAIW 1 


9 



4. How well does that distance represent the separation 


on the earth? 

"5. How well can the transmission of a photograph be 

measured? 


the earth? 


6, How is that transmission related to radiance at 


Clearly, any working definition of "image quality" adequate 
for the task at hand must be capable of obtaining quantitative answers to 
such questions and must relate those answers to measurable properties 
of elements of the image processing system. 


2. 3 System Element Performance Characterization and "Image Quality" 

This section delineates the techniques commonly employed 
to characterize equipment performance and image quality. The purpose 
is to establish some concepts that will subsequently be used and to 
point out why some others are unsatisfactory for the study objective. 

In order to ensure familiarity with the current state of 
the art, a literature search covering the period from 196 8 to the present 
was undertaken. 

Four overlapping categories; optical transfer function, 
resolution, noise, and subjective 'image quality, provide convenient 
areas for discussion. 


2.3.1 Optical Transfer Function 

The optical transfer function has been shown to be a 
useful tool to characterize the performance of many imaging devices. 
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The blur introduced by an optical system can be 
characterized by a "point spread function" s(x) defined implicitly by the 
convolution: 


i(x) = 0(X) * 5 (x) {2-1} 

where: l (x) - image brightness, o(x) - object brightness, and = 
position coordinates in plane orthogonal to optical axis, (* denotes 
convolution). 

Thus: 


I 0 v ) - 0 (v) z (is) 


( 2 - 2 ) 


where: l(v)^ i(x) , o(v) ^ o (x) , T(js)^S(^) , and 
denotes a Fourier transform. 

tT {'7) is the optical transfer function (OTF). It is in general a two- 
dimensional, complex valued function. \ z(v) j is the modulation 
transfer function (MTF). It is emphasized that the OTF is not a 

measure of image quality but merely the frequency response function 
of a linear device and consequently a measurable property of the 
performance of that device. 


Application of the OTF concept to photo -optical systems 
requires linearization of the generally non-linear development process. 
Generally the OTF concept is applied to the object to exposure image 

transfer process. Methods for measuring the OTF are discussed in 
Appendix A. 
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2, 3. 2 Resolution 

Resolution is related to the ability to determine object 
characteristics- especially shape from an 'image. A number of resolution 
criteria are in use. 

Rayleigh's criterion assumes that the diffraction limited 
images of two points are just resolved if the central maximum of one 
lies on the first minimum of the Airy disc of the other. 

The Rayleigh criterion is clearly related to the spread 
function and can thus be derived if the OTF is known. 

A number of resolution criteria are simply defined by 
an observer's ability to distinguish the existence of a particular target. 
Such criteria depend, not only on the properties of the image, but on the 
properties of the detection process as well. The most common is the 
standard Air Force Tri-Bar target. One can obtain a "modulation 
detectability curve" by having a number of subjects observed tri-bars 
of varying spacing and contrast and plotting the detection threshold 
contrasts versus spatial frequency.. "Resolution" is then defined as 
the intersection of an MTF and a modulation detectability curve. 
Uncontrolled variables and experience produce uncertainty in this 
measure of resolution. One summary measure of "image quality" 
that Is in use is the area enclosed between the modulation detectability 
curve and the MTF. 


The term "resolution" is sometimes applied to the ability 
of an optical system to "resolve" a specified object. This definition is 
similar to the preceding one but requires recognition as well as detection. 
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2.3.3 Noise 


In photo-optical systems, the major noise source is the 
granularity of the emulsion and is expressed in the granularity constant, 
G. 1 The rms density fluctuation observed in scanning a uniform density 
area is ^ 

•TaI 


where: a a = area of scanning aperture 


(2-3) 


In actuality the'emulsion records the continuous exposure distribution 
as a discrete, thin but nevertheless three-dimensional, distribution 
of silver particles. The photographic macro- image is a continuous 
intensity distribution which results from multiple scattering of 
photons traversing the developed emulsion. If the photographic image 
is observed over a very small area very close to the emulsion surface, 
it is not clear how the observed intensity is related to the intensity 
distribution which exposed the emulsion; in other words, the micro- 
image is not yet adequately understood. Considerable research has 

4-£ 

been performed to attempt to characterize photographic granularity. 

It has been represented as both additive and multiplicative noise. 

Since no clearly superior model exists It is most often represented 
by additive white gaussian noise. We employed this approach and 
represented the grain noise by rms fluctuations in transmission. 


For the electronic image processing system elements, 
an additive white gaussian noise model is theoretically as well as 
pragmatically acceptable. 


2. 3. 4 Subjective Image Quality 

Efforts have been made in a number of studies 
to define subjective assessment of image quality in a quantitative manner. 


ORIGINAL PAGEIS 
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Such techniques by definition include human variables which are not 
well controlled. It is not surprising that a universal subjective meas.ur.e 
has not been accepted although correlation of subject response with 
measurable parameters within the limits of specific product use has been 
shown. 


Subjective image quality efforts are directed to 
achieve a causal relationship between measurable system element 
performance properties, such as frequency response (OTF), signal-to- 
noise ratio, etc. , and the ability of the user to make subjective 
judgments (usually in the form of detection/recognition decisions) on the 
output product. The motivation of such efforts is cons istent with the 
objective of this study. However, the quality measures depend on the 
human detection process as well as the image data system. To develop 
quality control procedures the selected measure should depend only 
on the syctem itself. Thus, subjective image quality measures are 
not appropriate for the present study. 


2 . 3.5 Conclusion 

For most conceivable image processing systems, the 
elements' performance can be characterized by OTF's (if linear), 
nonlinear gains (photographic development), noise sources, or combinations 
of the three. Thus adequate "performance" descriptors exist. But 
based on the survey of the state-of-the-art of image quality evaluation, 
no objective technique to relate such descriptors to a quantitative 
measure of image quality was available. 
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SECTION 3 


MATHEMATICAL BASES 

The review of image system performance evaluation 
techniques indicated that descriptors of the influence of system elements 
on an image exist and are in common use. In particular, the performance 
of elements linear in intensity can be described by OTFs, the performance 
of elements non-linear in intensity can be described by non-linear gains, 
and elements which contribute noise can often be described by gaussian 
statistics. In general, the influence of a piece of hardware can be 
unambiguously described by a sequence of mathematical operators 
corresponding to these three measurable descriptors. 

One finds no such consensus on a metric for the quality 
of the output product of an image processing system. 

What the user of imagery does is to make observations on 

the output product, from which he estimates a radiant distribution (as a 

function of position) on the earth. But because the system is both band- 

limited and noisy, these estimates will necessarily be imperfect. A 

suitable image quality .metric is therefore the error which the system 

introduces into estimates of ground radiant distributions. Precedent for 

12 13 

this idea may be found in a study ’ of Lunar Orbiter imagery. 

In this section, it is shown that to formulate the image 
quality problem as one of parameter estimation yields an objective 
mathematical relationship between output product quality and element 
performance descriptors. In addition, analytic relationships required 
for implementation of the approach are derived. 
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3.1. Cramer -Ra'o Bound 


Consider a user who wishes to measure some parameter 
"A", say the radiance of a wheatfield, as a function of position "x" on the 
earth by making observations on an NDPF output product. Let the 
ground object be mapped onto the photograph space as the function S ( Xj . A), 
Since the system adds noise to the signal, the user will obtain 
a distribution of values with some standard deviation 0^ . To avoid 
consideration of the user's specific technique, it is assumed that he uses 
an optimal measurement process. 


It can be shown (see Van Trees “i that the variance, C^ 2 , 
of an unbiased estimate of A made from an output signal contaminated 
by additive white gaussian noise is bounded by the Cramer -Rao inequality: 

(3-1) 


2 

cr > 




-°\ 3SpL A) 
. A 


-.2. 


d x 


„ where N o = white noise power density (double-ended^) at the output and 
% g - record length. 


The user will generally be interested in estimating 

multiple parameters. Equation (3-1) can be generalized, for the 
"n" independent parameter case to: 




(3-2) 


* Van Trees uses a single-ended spectrum which results in a factor 
of 2 in the denominator of Equation (3-1). 
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where the matrix M is given by: 



dS(Z,A) 
d Ay,' 


A) 

3A J 


(3-3) 


A is simply the vector whose components are the independent parameters: 
A^, . . . , used to describe the input signal. 


If the noise power spectrum, is not white at the system 
output, an inverse or "pre-whitening" filter must be included as the 

last system element in order to apply the bound. 


The equality in the error bound represented by Equation 

(3-2): 

(1) holds if s(X,A) is linear in A 

_s- . A * 

(2) is approached when s(X } A) is nonlinear in A , but the 
signal-to- noise ratio becomes large. 


If one can compute the set of partial derivatives 


35(Z,A j 

3 A c 


as a function of changes in 5 (% , A ) resulting from changes in the 
performance descriptors of elements internal to the system, one can 
use this relationship to predict the change in estimation error intro- 
duced. Equation (3-2) is therefore the objective relationship between 
"system element performance" and "image quality" required. 


In order to illustrate the use of this tool, some specific 
examples are given: 


First, consider the task of. determining the boundary 
between two different crops and the reflectance level of each. The ground 
object can be represented as a three parameter "edge" target where A^ 
is the radiance level of one field, A^ the radiance level of the other, and 
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Ag as the location of the boundary. One computes the sensitivities: 

i -1, 2,3, determines the components of M from equation (3-3), 
3 A c — i . 

-obtains- the- inverse- -matrix- M ' , -extracts the' diagonal elements, and 

obtains the estimation errors on each parameter from equation (3-2), 


If the user utilized an optimal measurement technique, 
he would obtain statistically distributed values for each parameter with 
standard deviation <T£ . . according to (3-2). If his technique is not 
optimal, he will do worse, but in no event will he obtain more precise- 
results . 


As a second example, consider measurement of a river. 
The river can be represented as a four parameter target with the 
radiance level of the land on both sides, A^ the radiance of the river 
water, A^ the width of the river, and A^_ the location of the river. 

An example similar to "resolution” might be two bright, 
narrow lines on some background of radiance A , separated by distance 
A^. The "resolution" would then be expressed by er^ as a function of 
A ^ and A . 


It can be shown that a very broad class of functions 
(ground objects) can be expressed as n-parameter sets. Thus, 
this approach can, in principle, be applied to any conceivable ground 
object. It is sufficient for the purposes of this effort to use simple 
targets with four or less parameters shown in Figure 3-1 and discussed 
above. 
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Note that the estimation errors are directly 
proportional to the white noise spectral density, N q , the height of 
the double-edged noise power spectrum. Thus in the limit oi.a. 
nois"ele s s system, the object can be measured exactly regardless 

of the effect of the system. (It is of course implicit that the object 
form is a priori known in that it has been represented as a complete 
orthonormal set of n-parameters and that the system effects likewise 
are known a priori. ) 

In summary, a procedure has been defined for objective 
determination of the effect of changes in subsystem element performance 
on the precision to which users can make measurements of target 
characteristics from image products. 


3.2 Dimensional Analysis 

The dimensional analysis is useful in augmentation of 
intuitive understanding of equation (3-2) and of some of the results to 
be presented in Section 6. 

The one-dimensional case is considered. The ground 
object (parameterized target) is defined by radiance as a function of 
position. Call the radiance unit Jr| and the position (distance) unit 
|X| . The autocovariance function f' (x) is defined by 

U7n -P ] 

where f (x) in the image of the ground object and consequently f' (x^has 
dimension j R j 2 . 

The double ended power spectrum: 

, cCztttPz 1 ) 

H = j f «>e </*’ 

-oo 

2 

thus has dimension [ R | |xj . 


where j) = spatial frequency 



Therefore, 



, Now if the inpat target is expressed in terms of a set of 
orthogonal parameters A^, the the image is given by ?0c) •='5Cx, A ) 
for some particular A. By equation (3-2), the variances characterizing 
the error in estimating the components of that particular A vector are: 





. 



j 




By definition 


so that 



Finally, since 





substitution yields 



And the variances are seen to have the correct dimensions . 
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3,3 Noise Characterization 


A major noise _contribu.tp_r _in many-ima-ge-data -systems’ 
is photographic granularity. The photographic granularity constitutes 
a two-dimensional noise field which, for the purpose of this study, can 
be characterized by a two-dimensional white noise spectral density. 
Transformation relations" hips between the two-dimensional field and a 
one -dimensional representation are required for two reasons: 

1) Many users utilize microdensitometer traces or 
scans to obtain quantitative measurements from photographs. The 
microdensitometer record contains one-dimensional noise obtained 
by moving an aperture along some line over the two-dimensional 
photographic noise field. How do we account for the two-dimensional 
nature - of the image and include the appropriate noise spectral density 
in the calculation of the Cramer-Rao. bound (Eq, (3-2))? 

2) Section 4 will indicate practical reasons for a one- 
dimensional system simulation. Such' implementation requires valid 
characterization of the noise field in one dimension so that the 
simulation's results converge to the real physical situation. What 
magnitudes of rms noise should be added to the one-dimensional 
signals in simulating an image data system? 

The questions are related since the measuring aperture (size and shape) 
will determine the magnitude of the rms noise to be added in the simulation 
as well as the noise content of the one -dimensional trace used to obtain 
quantitative measurements (parameter estimates). There are two aperture 
shapes which might be employed: slits and circles. Both will be treated 

in the discussion that follows, however, a slit aperture appears to be the 
more likely candidate. 
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3.3,1 Treatment of noise for parametric estimation from one- 
dimensional scans. 

The procedure for establishing the variances of two- 
dimensional target parameters from one-dimensional scans requires 
knowledge of the scan noise spectrum. This, in turn, depends on the 
net two-dimensional spectrum as transmitted through the optical system 
to the point of scanning, and on the scanning aperture function. 


Arguments are presented illustrating that slits are a 
scanning aperture compatible with both the spirit of the present study 
and a mathematical requirement for the validity of the variance estimator. 
The results for a circular aperture are included for completeness. 

As either the slit length becomes arbitrarily large, or the width 
arbitrarily small, certain simplified asymptotic behavior between 
the scan spectrum and the two-dimensional field spectrum emerges. 

We can refer to such apertures as "partially asymptotic". In the 
event both conditions are met (totally asymptotic slit), the following 
relationship occurs: 


where: * 3T.e,(^s) 

1 


§ s (0 = -j §(A>°) 

- scan spectrum (transmission^-mm/cycle) 

2 2 2 

- field spectrum (transmission -mm /cycle .) 

- length of aperture (mm) 


')h 1 - spatial frequency in scan direction (cycles /mm) 
- spatial frequency orthogonal to 2^ (cycles /mm) 


* All spectra are double-ended (i. e. , - oq< f) <oo ) 
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The remainder of this section: 

• Derives a convolution theorem having general utility 
to the following sections. 

-• Establishes the relationship between film granularity 
and the two-dimensional Wiener spectrum. 

• Derives the relation between the scan Wiener spectrum 
and the field spectrum for slit and circular apertures and examines 
various limiting cases for the slit aperture. 


3 . 3 . 1. 1 A Convolution Theorem * 

Consider a two-dimensional stochastic variable T(x, y) 
representing point transmission of a field. Allow this fieid to be con- 
volved with an aperture function A to yield a new field having point trans- 
mission T^. We write: 

T a (5,W = JtCk + S, y+y) (3-5) 

— oo 

with the conditions: 

E{l) = 0 

oo 

J| AU,<y) = 1 

- OC7 

Since we shall be interested in the Wiener spectra of the two-dimensional 
fields, the autocovariance (eC 3l 6)J of is computed: 

?> s (xj) 5 E [t a (4 + X , Tj *P) r A ( 4 , 7)1 ( 3 - 6 ) 


*An equivalent theorem in frequency space is given by O’Neill in 
Reference 15. 


24 



which along with Eq. (3-5) results in: 




O 

| £ jVot + § , y +tj) T(x'-»-$ + eC, y'+y+fijj A(z,y) A(z\y'j d%dyd% dy 


or: 




°o 

frri' 


JO 


^(z'~x-k> c } y'~ y-^/3) A(x,y)A(z',^') d x dy, 


(3-7} 


-CO 


where ^ is the autocovariance of T. A change in variables % -*JJ. , 
according to: 


/c = £ - ; Z + cC 

allows Eq. (3-7) to be written: 


& = 


<£(/*., at ) 


-CO 


c*o 


d/x. d.Ar- ( 3 - 8 ) 


a(x, y ) A ( pL + X ~ cc , nr + y - fi ) d * dy 

) 

The integral in brackets is the auto-”covariance n * ( ) of the aperture 

and we having finally: 

<=o 

p s (oC } fi)= | f(pL-rK. % AT+0) d/1 4 nr* 


(3-9) 


This convolution theorem will have general utility to our analyses that follow. 


3,3,1. 2 Granularity and the Two-Dimensional Spectrum 

The NDPF model requires the Wiener Spectrum of photo 
graphic grain as inputs at several stages. 


This spectrum has already been assumed to be white 
thereby allowing its single degree of freedom to be evaluated in terms 
of the classical measure known as ’’granularity". Granularity is the 
RMS fluctuation in density ( cr^ ) derived from a scan of uniformly exposed 
film. For this purpose the scanning aperture is always large and the 


* A(x, y) is, of course, not a stochastic variable. 
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corresponding fluctuations are small. We shall deal here with the RMS 
fluctuation in transmission ( <X) T _ ). rather than -dens-i-t-y, but the- two are'' 
simply related. 


The variance of a one-dimensional scan is given by 

its autocovariance 0<,C°O (- (oc,cOj evaluated at zero. Thus we have 
through Eq.(3~9): 




= (o, o) 


<P A (/c,/lr) <LjxcLaT- 


'N 


Since the grain spectrum is white its autocovariance is written: 


whe r e : 


S O) 6 (nr) 




- two-dimensional Wiener 
{transmission - me ter 


spectrum of grain 
^ / cycle ) 


2 

and the variance 0"^. becomes: 

I 5 A (°»°) (3-10) 

A uniformly transmitting aperture of area A is written: 

A 


f l j A a (cucf/ltTt A A ) 

0 (outstcLe. A a) 

This yields 0 A (o,o} = iyA A and the value for in Eq. (3-10) is established 

Jf ■- A a (3-U) 

This is equivalent to Equation (2-3) presented earlier. 
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3. 3. 1.3 Wiener Spectra Relations 


Prior to deriving the relation between one and two- 
dimensional spectra, it is appropriate to digress in order to present 
the rational for selecting slits as having special significance to the 
present study. We build up the logic as follows: 

(1) The purpose of the present study is to determine the impact 
of processing parameter errors on optimum user perform- 
ance and not user performance per se. To this end simple 
psuedo-realistic one -dimensional targets (such as edges) 
are sufficient to make the connection since they contain 

the fundamental properties of transmission level and 
position. 

(2) In light of both (1) above and the point that the variance 
estimation represents optimum information extraction, 
the scanning aperture should not further degrade the 

•signal. Thus, its dimension in the scan direction should 
be small. 

(3) Generally, the lower bound of Equation (5-2) is 
realized only for differential departures of the parameters 
from their true values. The validity of this requires 
small noise fluctuations which in turn means a scanning 
aperture of large area. Since its dimension in the scan 
direction must be small, the other dimension must there- 
fore be large. Since the targets are one -dimensional, no 
signal degradation will occur under this requirements. 
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We now derive an expression between the scan spectrum 
resulting from a slit and the two-dimensional field spectrum. First, 
the scan spectrum <£ 5 (Vec) is defined as the Fourier transform of 

the. scan autocovariance <£> 3 (<x) : 


OO 


C^oc) ~ 


<t> s (oc , o) e 


-zirt^-cC 


6 cc 


- Oo 


Eliminating hy means of Eq, (3-9) yields; 


£0k) = 


t- 

“ C <3 


$(pL-t-oz , nr) (al , nr) d/JL <Lnr\ 
The field spectrum $ is defined by; 

<p O;*'') = 


-iTTtV^ 

e. c/cc 


OO 

n 


3? (}?/ JL^Vnr ) S cL V/tr 


(3-12) 


Eliminating ^ from Eq. (3-12) allows the integration over to be 


carried out in closed form: 

OO 


= 


— CO 


I e 






V 


^cXu-^nr) d/JL <Ljir ( 3 ~13) 


The aperture autocovariance, (p A , defined inEq. (3-8) is written for 
a slit of width pj~ (in the scan direction) and length 1; 



' ur 


W) = A (within A A ) 


(otherwise) 


Inserting in Eq. (3-13) we have: 


*» 0 W-ih 


MT 

t 7jut 







3 



i § C^ccAhr) 


% 

1 

OO 

* 

u 


ZjTtVtrV- 

fynrQv) e ^ 
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The two transforms in square brackets have solutions: 


I 

~MT- 

l 


\ 


t XI = si-nx. (rr jjj-) 


(X) - £ SC TIC '(7 T^l) 

U • J 

Thus, we have the final general relationship: 


J?5 VJ = <jr\> 0 c m^)\ $ S4% c CTTP^JL) dL^ 


(3-14) 

GENERAL 


We now give five special cases of Eq. (3-14) ? a n 0 f v/hich follow .immediately. 
The first three are the asymptotic slits defined earlier: 


$*<»«)* $QU>% -) S674C (JTXTltr) cLV a 


tr 


(3-15) 

NARROW SLIT 
( w ~^>o) 


£*(Hc) = jr sc^cCrr^Vcc) $ Ok , o) 

^ s ok)= i" ^ ol, °) 


(3-16) 

LONG SLIT 
(i oo) 

(3-17) 

LONG AND NARROW 
SLIT 


Finally, when both the dimensions of the slit become small or lar^e 

© 

together we have: 


OO 

=. [ i C Ac. ; T-L-) c/tL- 

- CO 


(3 -IS) 

SMALL SLIT 



S 6 7SC 


(3-19) 

LARGE SLIT 
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The asymptotic form given in Eq. (3-17) is recommended since It represents 
optimum parameter estimation. However, formally requiring the field 
spectrum $ at all frequencies can be side-stepped in those cases where 
the signal spectrum is severely band-limited. Unde.r .such -a- condition, 
the spectrum can be taken as white over the frequency domain of 
the signal: 

I Ok , °) — $ °) 


withy/ now interpreted as the net grain spectrum. Thus the scan noise 


spectrum is simply: 




(3-20) 


Now the field spectrum^V* is given directly by Eq. (3-11), and we have; 

o~r 


A *o = 


i 7 


(3-21) 


where (J). is the composite granularity due to all contributing system elements, 

An expression for the scan noise spectrum for a circular 

16 * 

aperture was developed during an earlier-study by Trabka , namely 


^5 fa) 


4£ ^(Srrd c v*) 

(rrd 0 Vo:) 2 


(3-22) 


where c/ c is the diameter of the scanning aperture and Trabka 1 s normalizing 
constant k = 4/c/ 0 . The aperture diameter must be large enough to vali- 
date the assumption of white noise made by Trabka. .Assuming that 5^ fa) 
varies slowly over the signal spectrum (effectively white noise) we set 

*s(°) -- N'-Wt £ 


since H 2 (Zl)/z 2 ~ 8/ 57r for 2*<1. 17 Using Eq. (3-11) for // we have 


*The units of the independent variable have been changed from radians /mm 
to cycles/mm in comparison to Trabka's Equation (6). 

30 



(3-23) 


_ 32. 

N o ~ 3*2 


<*r a a 




1,05 


CTr 

d 0 ' 


which is analogus to Eq, (-3-21) developed for slits. We note that both 
equations can also be expressed in terms of the rms density <T 0 through 
use of the approximation 

{T r c; 2.3 T 0 f 0 . (3-24) 


3.3.2 RMS Noise for One-Dimensional Simulation 

If the image data system is represented by one -dimensional 
model we must determine an equivalent rms noise, , to be added to 

the signal to represent the film granularity or noise. Since we simulate 
at intervals, , the simulated signal is band -limited by the Nyquist 
frequency, , we should only add noise inside this bandpass. 

In this case 


eg 2 = N 0 ( 2 ) = 


Ho 

A 


(3-25) 


We can relate to Kodak measured grainularlty by using Equations 
(3-21) or (3-23) depending upon the shape of the aperture we designate 
in scanning the image to obtain the one-dimensional trace. Consequently, 




2 _ 


* 1 



(long and 
narrow slits) 

{ circles ) 


(3 -26a) 
(3 -26b) 


If Kodak’s published granularity values are used for Cr , then A a in 
Equations (3-26a) or (3-26b) is the area of the aperture used by Kodak 
in making the measurements. 
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SECTION 4 


SOFTWARE IMPLEMENTATION 

To use Equation (3-2) requires calculation of the partial 
^ c / ^ A ) 

derivatives ___> I . There are doubtless image data systems 

d Ai 

for which the system can be represented analytically and the derivatives 
of the output signal obtained directly. But for a multi-element system 
which contains nonlinear elements and where it is desired to arbitrarily 
vary the response of individual elements, an analytic approach is 
intractable. Consequently, a computer program was written which 
simulates the system in order to calculate the derivatives and the 
Cramer-Rao variance bound. 

4. 1 Image Data System Simulation (IDSS) 

The primary function of the Image Data System 
Simulation (IDSS) is to compute the Cramer-Rao bound on the 
parameter estimation error introduced by a multi-element image 
data system composed of an arbitrary sequence of linear elements, 
nonlinear gains, and additive noise. A second capability of IDSS is to 
provide graphic display of test signals as they appear after every 
element of the image processing system so that the effect of individual 
elements can be visually assessed. A third is to determine the region 
of approximate linearity of the lumped-multielement system which 
contains individual nonlinear elements; the bounds of the "small signal" 
or "low contrast" linearizing approximation can be determined. 

Although IDSS was developed as a tool necessary for 
solution of the quality control problem, application is in no way 
restricted to the NDPF system. The program could equally well be 
used in analysis of entirely photographic, electro-optic, or digital 
systems. 


Program IDSS is written in Fortran IV for use on 
Cal span' s IBM 370/168 computer. It requires 250K core and of the 
twenty three subroutines, eleven are called from the 370 system 
library. 
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What the program actually does is create a numerical 
isomorphism to the image processing system being studied, create test 
signals, propogate those signals through the simulated system, and 
perform the mathematical operation on the processed output signal 
which yield the estimation error variances. 


The system simulation is one dimensional. A two 
dimensional simulation would involve no conceptual difficulties, but 
would greatly increase program size and running time. The one 
dimensional simulation is adequate for most conceivable purposes, 
but can given rise to certain subtle difficulties regarding the two 
dimensional photographic grain noise field which will be discussed 
in the second part of this section. 


Input 


• Most image data systems can be represented 

by a sequence of elements whose performance is characterized 
by either OTF's (linear elements), nonlinear gains, and/or additive 
gaussian noise. The data defining the response of each element 
is required as input. 


For the linear elements, the spread function rather than 
its Fourier transform, the OTF, is used for computational convenience. 

The spread function is of course always real valued, complex valued 
OTF's being manifested as asymmetrical spread functions. The effect 
of each linear element is thus determined by performing the 
convolution of its spread function with its input. The spread function 
defining the performance of each linear element is read as input. 

The nonlinear gains encountered are typically the 
photographic H-D curves. A table of output density versus input density 
values defining each nonlinear elements H-D curve must be read as 
input. The input signal is converted into density space (= -log (signal)), 
transformed to output density by linear interpolation of the table, and 
converted to transmission by exponentiation (10 -D 0ut ). 
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The noise is characterized by the standard deviation, 
which for the photographic granularity must be computed from Equation 
(3-26). A standard random number generating subroutine is used, its 
distribution scaled to the correct standard deviation, and the resulting 
white Gaussian noise added to the signal. 

Any input target type can be generated that is 
expressible as a function of position and a set of orthogonal 
parameters. The number of parameters and the magnitude of 
each is required as input. For- this study, targets with four or. 
less parameters are adequate. Consequently, the routine that 
creates targets, Subroutine GEN, currently produces: 1) One 

parameter target - constant level A(l), 2) Two parameter target - 

two delta functions separated by distance A (2) on background level A(l), 

3) Three parameter target - step at location A (3) between levels A{1) 
and A(2), or 4) Four parameter target - square pulse of width A{3) and 
location A (4) on background level A(l) with pulse level (height) A (2). 

These targets were illustrated previously in Figure 3-1. 

In all computations, the signal level is represented on 

a scale of 0 to 1 and will henceforth be referred to as ’'transmission’' . 

The 0 level corresponds to zero ground radiance. The 1 level corresponds 

to the radiance input which yields full- scale output voltage for the 

particular active sensor. ' This scale factor is defined as , Therefore 

all standard deviations in estimation of A^ representing radiometric 

parameters are in units of R m » The 0 to 1 "transmission" scale is 

actual photographic transmission at stages where photographs exist, 

« 

Output 

The program operates in any one of three modes. In 
all modes of operation the values of computational parameters and the 
functions defining the response of each simulated system element are 
-printed out. In mode 1 operation a table of the values and a table of • 
errors for each of the target parameters, A (I), is printed out. The 
process is repeated for new system element performance data as desired. 


34 



In mode 2 operation, the noise power spectral density is computed, and N 0 
and the pre -whitening filter spread function are printed out. Mode 3 
operation yields a plot of the target function as it appears after passing 
through each system element. Punched output of the processed target 
for auxiliary analyses (e.g. , edge gradient spectral analysis) is 
optionally available in any mode. 

Program Oper ation 

For logical convenience, the simulation separates the 
system into a number of blocks each containing sequentially a spread 
function, non-linear gain, and additive noise. A set of control cards 
determine whether one, two, or all of the three possible element 
types are present in any particular block so that any permutation of 
element types is achievable. A second set of control cards determines 
whether data is to be changed for each element so that the program 
user can vary the performance of any element or any combination of 
elements as desired. 


Mode 1 

A block diagram showing the major steps in the variance 
calculation mode is given in Figure 4-1. Control parameters, constants, 
and system element performance data are first read in. Each spread 
function is normalized to unit area; the spread function is integrated and 
then divided point by point by the value of the integral. Normalization 
yields conservation of energy (unity gain) in the subsequent convolution 
operations. Values of each of the A (I) are set and a test target generated. 
This target is propagated through the simulated system by the correct 
sequence of convolution with spread functions, and alteration by non-linear 
gains. Additive noise elements are omitted in this mode. The resultant 
output signal is stored. The operation is repeated for the same target, 
but with each of the A (I) in turn altered by a differential increment. The 
partial derivatives of the output signal with respect to each A (I) are then 
computed by subtraction of the first stored output signal from each of the 
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respective differentially altered output signals and division by the 
respective differential elements. The derivatives thus computed are 
labeled DYDA(I, X). The program next evaluates the integrals over 
X of DYDA (I, X) x DY DA ( J, X) for I, J=l, , . . , NA where NA is the number 
of A (I). The matrix thus created is then inverted and diagonal elements 
of the inverse matrix are extracted. These elements are labeled 
SIG(I, I). Each is multiplied by the white noise power density, yielding 
new values replacing the old SIG(I,I). The quantities SIG(I, I) thus 
arrived at are the variances in the respective A{I). 

The program then returns, assigns new values to the 
A (I) and repeats the computations for as many sets of A (I) values as 
desired. 


When all target A (I) values have been exhausted, 
the program returns, reads new system element performance data 
where' desired, and repeats the foregoing operations for the new 
simulated system. It stops when no new system element performance 
data remains to be read. 


Mode 2 

This mode is used to evaluate the' noise characteristics 
at the output of the system. The initial data input is the same as Mode 1. 
A one parameter (constant level) target is generated. The target is 
propagated through the simulated system, in this case with the noise 
elements included. No differentials or matrix elements are calculated. 
Rather the autocovariance function of the output signal is computed and 
Fourier transformed to yield the power spectrum and N Q , The operation 
is repeated several times (the number of replications is specified as 
input), and the results averaged. The prewhifcening filter, given by 
the square root of the reciprocal of the power spectrum, is computed 
and inverse Fourier transformed to yield the pre-whitening filter spread 
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function. This spread function must be used in Mode 1 operation as the 
last in the system simulation sequence if the power spectrum shows that 
the net noise is non-white. The value of N q can also be. .used a-s- input 
to Mode 1 for the double-ended white noise power density. 


' As in Mode 1, the operations are repeated for different 
target parameters and system element performance data as desired. 

Mode 3 

This mode is used to obtain plots of the target function 
at the output of various system elements. The initial data Input is the 
same as Mode 1. Targets are generated and propogated through the system 
as in Mode 1. The signal is plotted as it appears after each system element. 
No operations are performed on the output signal. 

As in the other two modes, the operations are repeated 
for different target parameters and system element performance data 
as desired. 


User descriptions, including a listing and flowchart of 
IDSS program, are given in Appendix B. 
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SECTION 5 


NDPF MODEL 

The initial step in applying IDSS to an image processing 
system is to specify its data flow and construct a model of the system. 
We will illustrate this procedure using the NDPF. The data flow in 
the NDPF Bulk Processing or Initial Image Generating Subsystem was 
discussed previously in Section 2.2. In this section, the NDPF will be 
broken down into elements whose performance can be characterized 
by OTF’s (frequency responses), additive gaussian noise, or non- 
linearities consisting of signal dependent gains. The results obtained 
from IDSS using this model are presented in the next section. 


5. 1 Element Performance Characterization 
NDPF Input 1 

The input is videotape records of telemetry from the 
satellite sensors. The telemetry and recording process are linear 
and characterized by frequency response functions (OTF's). White 
gaussian noise is added to the signal in transmission and recording. 
The RBV and MSS are linear devices characterized by OTF’s and 
white gaussian noise. 

Playback ' 

The playback videotape recorders are linear and 1 
characterized by OTF’s and white gaussian noise. The D/A 
conversion of MSS data and subsequent filtering is similarly 
characterized. 

EBR Electronics 

EBR Electronics are linear and characterized by an 
OTF and white additive gaussian noise. 
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EBR Pre emphasis Nonlinearity 


In order to obtain linearity between input radiance and 
archival transparency transmittance, and because of. processing con- 
straints on the EBR film, signal dependent gain is applied to the 
electronic input signal to the EBR, 


Photographic Processors 

The non-linear photographic process is defined by the 
Hurter -Driffield (H-D) curve which can be specified as a plot of output 
density as a function of input density 

(5-1) 

where D = density and T = transmittance. Since the linear elements 
of the system are linear in intensity (radiance), they are linear in T, 
so that the H-D curve represents a signal dependent gain. Henceforth, 
the relationships defining the photographic processor performance as 
well as that of the EBR preemphasis nonlinearity will be referred to as 
nonlinear gain. 

Contact Printers 

The physical process involved in contact printing is 

coupled multiple scattering within the master -duplicating film sandwich. 

That process yields a nonlinear image transfer in the general case, but 

has not been adequately studied to be useful in the present application. 

As part of a study of the Kodak Niagara Printer, or.e of the authors (EKS) 

showed that in the limit of infintessimally thin continuous emulsions, .the 

contact print would be the Fresnel diffraction pattern of the master. 

18 

Fresnel fringes for edge targets were in fact observed by Yeadon. 
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Diffraction phenomena result from discontinuities in the 
electric field. However, as edge gradients become finite and sufficiently 
low, the diffraction effects become less pronounced and the printing 
process approaches linearity. For the sensor imposed bandlimit 
of approximately 50 cycles /mm sufficient edge blur is expected that 
the linear model be valid. Should future payloads have increased OTE 
cutoff frequencies, this contact printer approximation must be re-examined. 

In the present application, therefore, the contact printing 
process is considered to be a linear process with white gaussian noise 
(film granularity). The film MTF is lumped together with the printer > 

OTF and referred to as a film-printer OTF. 

Bulk Enlarger 

The Bulk Enlarger is a linear device characterized by 
an OTF. That OTF is expected to vary as a function of field position, 
but since its magnitude at the sensor cutoff frequency {^12 cycles/ 
mm on 9. 5 in. format) should be high, field position variation can 
be neglected. Barring catastrophic occurrences, the enlarger's 
aberrations are constant, and the source of OTF variation expected 
in normal operation would be defocus. (Note that if the enlarger 
system is misaligned, there will be defocus as a function of field 
position, which cannot be neglected. ) 

The film MTF. is lumped with the enlarger MTF so that 
the subsystem is modeled as a linear element characterized by its 
OTF, and additive white gaussian noise. 

5. 2 Model Constraints 

Several constraints were imposed on the analysis of 

the NDPF: 
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1. As the RBV sensor in the current payload has been 
inoperative, the MSS configuration of the image processing system is 
emphasized. 

2. The bulk processing (Initial Image Generating 
Subsystem) only is analyzed. Precision processing (Scene Correction 
Subsystem), the prime purpose of which is large scale geometric 
modification for mapping purposes, is not included in model manipulation. 

3. Paper prints are not likely to be used quantitatively 
and are not considered. 

4. The bulk processing is modeled only through the 
third generation transparency which is the prime product delivered 
to users for quantitative study. 

5. Sensor -telemetry noise can be lumped together, 
and NDPF electronic element noise can be neglected. 

6. Telemetry and recording apparatus can be assumed 
transparent (no effect on signal). 


5. 3 MSS Configuration Model 

In terms of the subsystem element performance 
descriptors discussed, and under the preceding constraints, the 


NDPF model for processing MSS 

Element 

1. Sensor 
Electronics Noise 

2. EBR Control Electronics 

3. Digital to Analog (D/A) 

conversion 

4. Bessel Filter 

5. EBR Electronics 

EBR Preemphasis non- 
linearity 


imagery is as follows: 

Measurable Performance Descriptor 
OTF 

Standard deviation (gaussian, white) 
OTF 

OTF 

OTF 

OTF 

Nonlinear gain 
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6. EBR (Beam) 

7. EBR (Film) 

EBR film development _ 
EBR film granularity 

8. Contact printing /enlarging 
Contact printing/ enlarging 

film development 

Contact printing /enlarging 
film granularity 

9. Contact Printing 
Contact Printing film 

development 
Contact Printing film 
granularity 


OTF" 

OTF 

Nonlinear gain (H-D) 

Standard deviation (gaussian, white) 
OTF 

Nonlinear gain (H-D) 

Standard deviation (gaussian, white) 
OTF 

Nonlinear gain 

Standard deviation (gaussian, white) 


This model is represented by a series of 9 blocks identified above. 
Input values for each of the performance descriptors are presented 
in Section 6. 1. 


5. 4 RBV Configuration Model 
Element 

1. Sensor 
Electronics noise 

2. EBR Control Electronics 

3. EBR Electronics 

EBR Preemphasis nonlinearity 
4-7. Same as MSS elements 6-9 


Measurable Performance Descriptor 
OTF 

Standard deviation (gaussian, white) 

OTF 

OTF 

Nonlinear gain 


5 . 5 Large Scale Errors 

If all of the system elements introduced no degradation 
into the signal, the ERTS/NDPF system would only introduce a scale 


OUTSWAP ? AGE ® 
OF POOR QPALirC 
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change into the ERTS imagery. However, effects of payload attitude, 
sensor optics, and distortion internal to the NDPF result in a film 
plane geometry distorted compared to the g-rou-nd. The ND'PF was 
designed to correct for this mapping distortion (correction .mechanisms 
are outlined in Appendix A. ); however, some residual geometric 
mapping distortion error results. This error applies to measurement 
of distances of the order of magnitude of the frame size. 

If the NDPF is described by a pseudo spread function of 
the .entire lumped data processing system, the geometric estimation 
error introduced by the NDPF element performance descriptors is 
separable from the large scale geometric mapping distortion described 
above; the two error sources can be treated independently. The errors 
calculated assess the small scale, not the large scale geometric error. 

Similarly the radiant input to the sensor would be mapped 
into a related radiance scale at the output product for an ideal system. 
However, a radiance mapping error exists over distances of the order 
of the frame size due to sensor and processing nonuniformities. The 
correction mechanisms for this error are also discussed in Appendix A, 
and the residual errors resulting from application of the algorithm can 
be calculated directly. Again, this radiometric error is separable from 
the parameter estimation errors of concern in this study. 

5. 6 Modeling Summary 

The NDPF and sensors have been modeled in terms of 
accepted, measurable descriptors of the performance of system elements. 
These descriptors can be utilized as a sequence of operators to obtain 
system output as a function of input to the sensor. To develop or assess 
quality control tolerances requires determination of the effect of changes 
in the performance descriptors on parameter estimation errors. 


44 



SECTION 6 
RESULTS 


The results of application of the technique to the NDPF are 
presented in this section. The first part discusses the input data base 
identifying the specific values assigned to the descriptors (functions) 
which characterize the performance of each element of the NDPF model 
delineated previously. Tn the subsequent discussion, the nominal 
performance of the NDPF/MSS configuration is established and 
examined in detail. The effects of off-nominal element performance 
are determined. The results expected from the RBV configuration 
are considered by comparison to the MSS, The implications of these 
results on procedures for quality control of the NDPF elements are 
subsequently examined. Finally, a technique for calibrating ERTS imagery 
in order to relate transmittance to ground radiance or reflectance is 
illustrated. 

6. 1 input Data for Element Performance 

Valid quantitative definition of each spread function, 
nonlinear gain, and noise standard deviation are obviously required 
as input to IDSS. Approximate accuracy of input data is adequate 
to determine the influence of changes in system element performance 
or changes in estimation error, which is of course, the prime objective. 
However, accuracy of the calculated estimation errors as a measure 
of absolute system performance requires high confidence level input data. 

Although some simple measurements of element 
performance were made by us, for the most part it was necessary 
to rely on external information sources. Sufficient confidence in 
the approximate accuracy stipuated is asserted; the best data 
found to be available was used. But.it must be emphasized that 
the reader keep in mind that inferences of absolute system 
performance are affected by absolute accuracy of input data. 

«8£*3iS£ 




Input data used and the manner in which it was 
acquired for each system element are given below. All spread 
functions plotted have been normalized to unit area. The unit of 
length used in computations is the micrometer and distances refer 
to the 70 mm format. All nonlinear gains are presented as plots 
of input density vs output density. Noise sources, always gaussian 
and white, are characterized by the applicable standard deviation. 

a) Payload Sensors 

1. MSS - the sensor spread function was obtained 

2 

by inverse Fourier transform of the MTF given in the User's Handbook, 
No phase data were available. The curve given appears to have been 
fitted to only four data points. The resulting MSS spread function is 
plotted as Figure 6-l(a-}. 


2. RBV - the sensor spread function was obtained by 

inverse Fourier transform of an "eyeball average" of the three 

2 

characteristic MTF's given in the User's Handbook. No phase data 
were available. The spread function is plotted as Figure 6 — 1(b). 

b) Telemetry and Recording - Assumed Transparent.^ 

c) Playback Video Tape Recorder - Assumed Transparent."^ 

d) Video Tape Recorder Control Electronics - Assumed Transparent.^ 

e ) EBR Control Electronics (EBR CTL) - Spread function obtained by inverse 

19 L 

Fourier transform of MTF given in ERTM-H-81. Plotted as Figure d~1(c) 

f) Digital to.Analog Converter (D/A) - As in (e) above. See Figure 6-l{d). 
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Figure 6-1 SPREAD FUNCTION INPUTS (PART 1) 
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■ g) Bessel Filter - As in (e) above. See Figure 6-l(e). 

h) EBB Electronics - As in (e) above. See Figure 6 -1(f). 

i) EBB Pre emphasis Nonlinearity - 

Two sourc.e.s. of. info rmation- were- used, 

2 . 

1) User's Handbook. The system is designed to obtain a net o = -1 
between log of the ground radiance and archival density. Further, 
the archival film is processed to T = 2 (with "slight 11 deviation from 
linearity at end of the H-D curve, the actual shape not being given). 
Thus, the EBB preemphasis nonlinearity has Y= Mapping the 
0-1 (relative) earth radiance scale onto the Handbook prescribed 
archival density range yields the curve shown in Figure 6-2(a), 

2) Measurements - A current quality control procedure utilizes a 
twenty-one step gray scale at the head and tail of each roll of film 
processed. A sample was obtained and measured, and an archival 
processor H-D curve obtained. Each frame of imagery includes a 
fif.een step gray scale written by the EBB for radiometric calibration. 
The densities were read on an archival frame. The relative log 
radiance of the gray scale was obtained from the. User's Handbook 
and an H-D curve for relative log ground radiance to archival 
density thus obtained. This curve together with the curve for 
archival processing alone allowed determination of the EBB 
preemphasis nonlinear gain plotted in Figure 6~2b, 

j) EBB Beam 

The EBB beam spread function was obtained from EBTM-H-8 1,*^ 

It is Gaussian in shape and plotted in Figure 6 -3a. 

k) Film spread functions 

The film spread functions were measured frc-.n sample data. The 
edge of the large square in a standard Air Force tribar target 
currently used in iquality control at the NDPF was scanned with a 
microdensitometer on archival, second, and third generation 70 mm 
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o .2 .4 .6 .8 1.0 1.2 1.4- 1.6 1.8 2.0 

DIN 

<b) EBR PREEMPHASIS NONLINEARITY (MEASURED) 


Figure 6-2 EBR PRE-EMRHASIS NON-UNEAR GAIN FUNCTIONS 
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(a) EBR BEAM SPREAD FUNCTION 



(b) FILM SPREAD FUNCTION 
Figure 6-3 SPREAD FUNCTION INPUTS {PART 2) 


50 


transparencies. The edge spectra were computed by an edge 

gradient technique described in Appendix A. Each spectrum 

was divided by that of the previous generation and the resultant 

MTF Fourier transformed to yield the spread function of a 

single duplication process. The MTF's were compared to the 

2 

limited section of the film MTF given in the User’s Handbook 
(Figure H-5) and found to be in agreement. It was concluded 
that the combined contact printer plus film MTF is approximately 
equal to the film MTF alone. Consequently the resulting spread 
function was used for both the archival film and for the second and 
third generation contact printer /film combinations. This spread 
function is plotted as Figure 6 -3b. Edge ringing (Fresnel diffraction) 
was not observed on the samples examined which provides support 
for the validity of the linear contact printer model approximation, 

1) Film Nonlinearities 

Three approaches were used: 

2 

1. User's Handbook. The design goal or nominal 
NDPF performance were obtained from the Handbook (Figure H-4) 

and are plotted in Figures 6-4 a-c for the first through third generations. 

2. The twenty-one spot quality control gray scale was 
measured on samples of each generation. The H-D curves thus obtained 
are plotted in Figure 6-5 a-c for the first through third generations, 

3. The fifteen level. calibration gray scale was measured 
for each generation. The H-D curves for the second and third generation 
processing thus obtained are plotted as figures 6-6 a } b. This 
approach could not be applied to the first generation as no independent 
measure of the EBR preemphasis nonlinearity could be obtained and the 
15 step density input to the first generation processor is unknown. 

The obvious disagreement between the three sets of data -will be examined 
later. 
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0 .4 .8 1.2 1.6 2.0 2.4 

D IN 

<c) THIRD GENERATION 

Figure 6-4 NOMINAL FILM H-D CURVES (HANDBOOK) 
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(a) FIRST GENERATION 



(b) SECOND GENERATION 



SECOND GENERATION DENSITY 
(c) THIRD GENERATION 
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Figure 6-5 MEASURED FILM H-D CURVES (21 STEP GREY SCALE) 
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(b) THIRD GENERATION 

Figure 6-6 MEASURED FILM H-D CURVES (15 STEP CALIBRATION GREY SCALE) 
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m) Electronic Noise 


Total electronic noise at the NDPF is assumed to yield 
electronic noise standard deviation equal to 10 . Electronic noise 

added by NDPF elements is neglected. ^ This noise amplitude is small 
compared to the photographic grain noise. 

n) Film Grain Noise 

The rms granularity values given in the User's 

‘ 2 

Handbook (Figure H~6) can be used to calculate £g , the rms noise 
fluctuation. The value of cTj A ^ , presented in Figure 6-7 for 
either the EBR (archival) or copy film, is substituted into equation 
(3-26) to compute . It is recognized that some noise sources encoun- 
tered may be signal-dependent; nevertheless, -t-he perturbation technique 
can still yield correct results if a standard deviation corresponding to a 
nominal or worst case is utilized. 

6. 2 Nominal Performance - MSS Configuration 

Referring to the MSS image processing configuration 
model established in Section 5. 3, the operations entered in a sequence 
of 9 blocks in IDSS to simulate the system are: 

1. Convolution with the payload spread function. 

2. Convolution with the EBR Control Electronics spread function. 

3. Convolution with the D/A Conversion spread function. 

4. Convolution with the Bessel Filter spread function. 

5. Convolution with the EBR Electronics spread function. 

6. Level change by the EBR preemphasis nonlinearity. 

7. Convolution with the EBR Beam spread function. 

8. Convolution with the film spread function. 

9. Level change by the nonlinear first generation development proces 
(H-D curve) 

10. Convolution with the film-contact printer spread function. 
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11. Level change by the nonlinear second generation development 
process (H-D curve) 

'12. Convolution with the film-contact printer spread function. 


13.- Level change by the nonlinear third generation development process 
' (H-D curve) 


Operation No. 10 is replaced by convolution with the Bulk Enlarger - 
spread function for 9. 5 inch products. Distances refer to the 70 mm . 
format unless explicitly stated otherwise. 

Recalling the disparity between the User's Handbook 

and measured H-D curve data, a decision was required on what to 

accept as "nominal". Since the measured data was taken from a 

small sample of NDPF products, and there is no evidence confirming 

the data as nominal or even typical, the Handbook data were chosen. 

This data is, however, clearly idealized, so the measured data were 

20 

also used as input ‘and the results compared. Other sources indicate 
that the first generation nonlinear gain or rr i is about 1. 5 in agreement 
with the measured data. 


In order to estimate the noise statistics the IDSS program 
was run in Mode 2 with nominal element performance data detailed 
in Section 6. 1. The noise spectrum was found to be white and a pre- 
whitening filter spread function need not be included in subsequent 
computation. 

To understand the noise propogation behavior in more 
detail, consider the last several steps in the data flow of the NDPF: 


© 

e 

© 

e 

« 

© 


additive noise (archival film grain) 
film -printer film spread function 
non-linear gain 

additive noise (copy film grain) 
film -printer spread function 
non-linear gain 

additive noise (copy film grain) 


L 2nd Generation 
Negative 


> 3rd Generation 
Positive 


■■ i t k ni?. 
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The input to these steps is a bandlimited, low-noise signal. Since 
the noise amplitude is small compared to the signal, the nonlinear gains 
do not modify the noise s.p.e.ctr-um. Now the film -printer OTF cutoff 
frequencies are large compared to the cutoff of the signal imposed by 
the sensor. The first grain noise is filtered by two OTF's, the second 
by one, and the third not at all. Small changes in the performance 
of either film or printer, will consequently have negligible effect on N , 
Therefore, for the current ERTS system, N q is independent of reasonable 
changes in NDPF eleme'nt performance. 

Another benefit is obtained from the insensitivity 
of N to NDPF system performance. From equation (3-2) 

°a\ * 

with the equality holding for the sufficiently high signal-to-noise ratio. 

The constancy of for slightly off nominal element performance allows 
relative error comparisons, namely: 

Cf? (off nominal) (M ' ).. (off nominal) 
t _ && 

p ____ i 

( X '• (nominal) (M' )^.. (nominal) 

2 _ 2 

where CTf - CT^ to simplify notation. Thus the quality control 

problem for the current ERTS image processing system can be studied 

by computing only the (M ^ )^. as a function of element performance. 

Therefore initially only the relative performance was examined. The 

effect of N was included later to relate relative performance to 
o r 

absolute values. 

6.2,1 Edge Target Results 

The three parameter "edge" target, where A^ = 
one radiance level, A^ = other radiance level, and A^ = the position of 
the edge between A ^ and A^ will be used to describe the change in image 
quality as a function of changes in subsystem element performance for 
the following reasons: 
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1. It is the simplest target that yields both radiometric 
and geometric estimation errors, 

2. It corresponds to a most common user estimation 


task. 


3. The spectrum of the output edge gradient can be 
used to determine system radiometric linearity. 


6, 2, 1. 1 Linearity 

In order to study system linearity IDSS was used to 
generate the nominal MSS output for the twenty targets shown in Table 6-1. 

The NDPF system output for each of these targets was 
subsequently used as input to an edge gradient spectrum computation 
program (Appendix D) and the spatial frequency response of the system 
thus obtained. 

If a system is linear, then its frequency response is inde- 
pendent of input. A nonlinear system on the other hand is characterized 
by harmonic generation, the presence of signal at higher frequencies is 
observed, and the apparent frequency response is input dependent. The 
region over, which the response of a nonlinear system is independent of 
input is defined as the "small signal" or quasi-linear region which is 
useful to those who want to approximate the nonlinear system by a linear 
one. 

The nominal NDPF is designed such that the four cascaded 
nonlinearities give a net processing Y = -1 so that ground radiance 
should be directly proportional to third generation transmission. How- 
ever, each nonlinearity generates harmonics while the O'TF of the followin 
linear element can eliminate some of those harmonics. Applying an 
inverse nonlinear function to the modified signal might not yield the 
initial input again. Therefore, the presence of the linear filters (OTF’s) 
between the nonlinearities may cause the system to behave in a nonlinear 
manner even though the net Y is held equal to minus one. 
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TABLE 6-1. EDGE TARGET PARAMETERS 


(100/t4.m records) 


Target Designation 



(normalized 
radiance 
. level) 



(normalized (edge location) 
radiance 
level) 


1 

« v 

2 

M 

3 

M 

4 

11 

5 

if 

6 

.25 

7 

u 

8 

ii 

9 

ti 

10 

ft 

11 ■ 

.45 

12 

If 

13 

u 

14 

1! 

15 

Jt 

16 

. 65 

17 

* ri 

18 

it 

19 

u 

20 

tt 


.0001* ■ 50. 

.2 

.4 

.6 

.8 

.0001 " 

. 2 
.4 
.6 

.8 

.0001 " 

.2 " 

.4 

.6 

.8 

w 

. 0001 " 

.2 ' ” 

.4 

.6 

.8 


(*The value . 0001 was used instead of zero because the program 1 must 


compute 



in the nonlinear gain subroutine.) 
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The edge gradient spectra computed are shown in 
Figure 6-8. Note that the spectrum of all the edges except those 
with the very low (. 0001 or . 05) radiance levels are identical and 
that the others all contain higher modulations as frequency increases. 

The system thus seems to be behaving in a linear manner over a wide 
but bounded target radiance range. From the results obtained this 
range is i at least, 0.2 to 0.8 (normalized radiance). 

In order to understand the cause of this result, the 
program was run in the signal tracing mode which yielded the plots 
shown in pages 214 to 233 of the data, Appendix C, corresponding to 
edges 1-20 respectively. Distortion of edges 1, 6, 11, 16 is obvious. 
Distortion can be seen in edges 2 through 5 as well, although it is 
less pronounced. From the plots, the distortion occurs at either 
the electronic preemphasis or first generation development non- 
linearities. Inspection of the EBR preemphasis nonlinearity curve. 
Figure 6-2 shows that the linear portion of the curve extends only 
as high as input density = 1.146 which is equivalent to a 
normalized input signal (transmission) of 0.07. Thus, the edges 
for which the system exhibits nonlinear response have at least 
one radiance level on the portion of the electronic preemphasis nonlinear 
gain curve which approaches the asymptote at (D.^, D^^) = (eo, .85). 

(In computation, infinity was taken as D = 1000). 

As the subsequent nonlinearity H-D curves are chosen 
such that any input signal which lies on the linear portion of the 
electronic pre emphasis nonlinear gain curve must also lie on the 
linear portion of all succeeding H-D curves, it is concluded that 
the system behaves in a linear manner if the signal' is not clipped 
by the 0. 07 breakpoint of the electronic preemphasis nonlinearity. 

If the four H-D curves do not indeed distort spectra, 
the spectrum representing the linear input range in Figure 6-8 should 
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Figure 6-8 EDGE GRADIENT SPECTRA FOR NOMINAL PERFORMANCE 


be matched by that obtained by the system with the nonlinearities 
removed (composite system OTF), 

The four H-D curves were removed, from the simulation 
model and the edges run again. The edge gradient spectra were computed 
as before and are plotted in Figure 6-9. The spectra for all edges are 
identical, equal to those in Figure 6-8, corresponding to edges 7-10, 
12-15, 17-20. 


It follows that the ERTS nominal system can be 
characterized by an end-to-end frequency response. One useful 
result of this conclusion is that if the NDPF is constrained to 
nominal performance, the OTF of the sensor in the payload can 
be monitored by measuring, on the output products, the edge 
gradient of a ground object. However, it is fallacious to ignore 
the nonlinear elements in computation of estimation errors as the 
following section shall show. 


Before proceeding, it is instructive to compare 
the edge gradient spectra of the input to that at later stages within 
the NDPF. The edge targets in the quasi-linear region were 
generated: 1) after the MSS, 2) after D/'A conversion, 3) on the 
archival product, and 4) on the third generation product. Edge 
gradient spectra were computed and are plotted in Figure 6-10. 

The D/A conversion is logically lumped as an integral part of the 
payload, even though it occurs on the ground. The NDPF processing 
is shown co modify the spectrum very little thus supporting the 
approximations made earlier in the report that were based on the 
signal bandlimit occurring before the addition of significant noise 
sources. 
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modulation 



SPATIAL FREQUENCY (CYCLES/mm) 


Figure 6-9 EDGE GRADIENT SPECTRA - COMPOSITE SYSTEM OTF 
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Figure 6-10 COMPARISON OF EDGE GRADIENT SPECTRA AT VARIOUS SYSTEM POINTS 


6 . 2 . 1.2 


Nominal Product Quality 


The nominal error bounds for- the- estimation of 
the three parameter edge targets were calculated as measures of 
product quality. The targets outside the linear range defined 
in the preceding section were eliminated resulting in a new 
set of 12 targets having a record length = 100 micrometers, a 
sampling increment = 0.5 micrometers, and a distance unit = 
micrometer. The target set is defined in Table 6-2. 


TABLE 6-2. REDUCED EDGE TARGET SET 


Designation No. 


1 

2 

j 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


(normalized {normalized 
radiance level) radiance level) 



(edge 

location) 


.25 .2 • 50. 

" .4 11 

" .. 6 " 

" .8 ” 

.45 .2 

" .4 ,r 

" .6 

" .8 n 

.65 .2 

" .4 " 

" .6 " 

" .8 " 
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These targets can be grouped according to the normalized radiance 
level difference V*/ ~ I a measure of target contrast. 

This results in the following classification: 



♦ 05 

♦ 15 
,25 
.35 
.45 
.55 


Designation No. 

1 , 6 , 11 

2, 7, 12 
5, 10 

3, 8 

9 

4 


IDSS computed the matrix elements (M 1 = 1, 3. (M ^ ^ and 

(M ^22' proportional to the error in estimating radiance levels, 
remained constant for all edges. They can be made arbitrarily 
small by increasing the record length- -it is the change in their values 
as a function of subsystem element performance for -a fixed record 
length that must be examined to evaluate quality control procedures. 

the edge location estimation component, was found to be 
independent of the absolute values of and A^, hut to depend on 
Z±R/R. m only. The results are plotted in Figure 6-11. The shape of 
the curve suggests the log-log plot of Figure 6-12. Note that Figure 6-12 
shows the data to be exactly colinear with slope of exactly -2. Now: 



so that the data demand: 


by equation (3-2) 


or 



where: ^ means T 'is directly 
proportional to" 
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This relationship is expected for a linear system, and can, in fact, 
easily be derived 'from equation (3-2) analytically. It confirms that 
this re.gion_.of input -target radiance levels represents a quasi-linear 
performance region for the NDPF. 


To determine the effect of the non-linear g_.ins, these 
steps were removed from the simulated system and the set of test 
targets was run again. The results are plotted In Figure 6-13 and 
6-14. Notice that the proportionality 


o' ~ 


A?? 


is preserved, but that the proportionality constant is roughly an order 
of magnitude lower. Thus we can conclude that -die presence of the 
four nonlinear gains increases the estimation error, even though the 
frequency response is not distorted. 


6. 2. 2 Square Bar Target Results 

Although edge targets will be used in examining off 
nominal system performance, the square bar target was also 
considered during the nominal performance analysis as an additional 
-descriptor. The following set of targets were generated where A^ = 
background radiance level, A-, = bar radiance level, A^ = bar width 
and A^ = bar location: 

Aj was held constant. A^ = 0. 1 

A^ was varied in steps of . 2 between A^ = .15 and A^ = . 95 
A^ was varied between A^ - 15 and A^ = 105 micrometers. 

A^_ was held constant at 75 micrometers, the record midpoint. 

Signal trace plots for these targets were shown in Appendix C, Pages 

234 through 278. Plots of the matrix elements (M *) , (M~*) , and 

~ i 2233 

(M ) as a function of bar width for each A^ value are shown in 

Figures 6-15 through 6-17. Both the radiometric and mensuration 
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WIDTH (MICROMETERS) 

Figure 6-16 VARIATION OF BAR TARGET WIDTH ESTIMATION ERROR WITH 
TARGET WIDTH AND CONTRAST 
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(M1CROMETERS/R 



WIDTH (MICROMETERS) 



Figure 6-17 VARIATION OF BAR TARGET LOCATION ESTIMATION ERROR WITH 
TARGET WIDTH AND CONTRAST 


S 


accuracies represented by these relative errors begin to increase 
for targets less than 40 micrometers and become asympototic for 
ta.rg.ets smaller -than -25 micrometers,' This corresponds to ground 
distances of approximately 134 and 84 meters respectively. From 
Figure 6-15 we see that the variance in the radiometric measurement 
is independent of target contrast, AR/Rm, while the variance in 
measurement of bar width (Figure 6- 1 6) and location (Figure 6-17)' both 
depend upon this contrast. If these data are presented on a plot of 
log (M versus log (AR/Rm) we find that both terms have the inverse 
square dependence on contrast described earlier in connection with the 
variance in locating edges. We note that the variance in target width 
measurement is an order of magnitude larger than the variance of the 
location measurement and is about twice the variance in edge location 
(Figure 6-11). The later result indicates that target width measurement 
is equivalent to two independent edge location measurements for sufficiently 
wide targets. 

6. 2. 3 Delta Function Target Results 

In order to obtain an estimate of nominal system 
‘'resolution", two parameter targets consisting of two narrow bars 
on a constant background were used. Signal trace plots of these 
results are shown in Appendix C, pages 279 through 283. The two 
bars are no longer distinguishable at a separation of 25 micrometers 
which corresponds to a distance of about 84 meters on the ground 
in agreement with the results from the square bar targets presented 
above. 


6,3 Effect of Off Nominal NDPF Element Performance 

The effect of variation in performance of the NDPF 
subsystem elements is now considered for the nominal (Handbook 
nonlinearity data) NDPF system. 
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'6. 3. 1 


Electronics 


The spread functions for the EBR Control Electronics and 
for the EBR electronics are practically delta-functions They are unlikely 
to vary and a great change would be required before the si'gnal was 
affected. These two spread functions were therefore not varied in the 
sensitivity analysis. 

The D/A conversion spread function is defined by the 
type of device used. Small variations in its shape are not likely to 
occurjeither it works or it doesn't. The same holds for the Bessel 
filter. Variation in the performance of these two devices is therefore 
not examined here. However, both devices are related to the payload 
performance --If a new, higher performance sensor were used in the 
future, both devices would also be redesigned. In this sense, they can 
be considered part of the payload. The effect of higher payload 
performance on product quality will be examined later. 


6 .3 . 2 EBR Beam 

The EBR beam can change shape due to variations in the 
electron optics, misalignment, contamination, etc. The effect of EBR 
beam width was examined. The set of test edges (see page 66) was run 
for EBR beam spread functions 20%, 50% and 100% wider than nominal 
shown in Figure 6 -3a. The off nominal EBR spread functions are shown 
in Figure 6-18* 


The effect of EBR beam width increase on Cf t and 

2 . 

is shown in Figure 6-19a. The effect on o' is shown in Figure 6-19 b. 
The graphs represent the data for all twelve test edges; the variation 
among the targets is shown by the vertical bar. 
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I i 

20% OFF NOMINAL 
(a) 




100% OFF NOMINAL 
(c) 


Figure 6-18 OFF NOMINAL EBR BEAM SPREAD FUNCTIONS 
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FRACTIONAL EDGE LOCATION VARIANCE ! , FRACTIONAL RADIANCE LEVEL VARIANCE 



PERCENT SPREAD FUNCTION WIDTH INCREASE 

(a) 



PERCENT SPREAD FUNCTION WIDTH INCREASE 

(b) 

Figure 6-19 EFFECT OF E8R BEAM SPREAD FUNCTION INCREASE ON LOCATION 
AND. RADIOMETRIC ERROR 

SSSSSSS 


2 

The variation in cr 3 is the largest. However, 
the standard deviation in edge location estimation for a T00% 

EBR beam width increase is only ten percent larger than the 
nominal value. This insensitivity indicates that abberrations causing 
assymetry of the beams and also minor beam misalignment need 
not be monitored. . The monitoring procedure need only be capable 
of detecting gross changes in the beam width or alternatively in the 
EBR OTF. 

6. 3. 3 Film-Printer Spread Functions 

• The sensitivity of the estimation errors to changes 
in the film-printer spread functions was next examined. The edge 
target set was propagated through -the system for 20%, 50% and 
100% increase in each of the three film stages (archival film, 
second or third generation printer-film). Figure 6-20 shows plots 
of the three off nominal film or film-printer spread functions. The 
effects on the estimation errors were identical regardless of which 
of the three elements were degraded. The results for all twelve 
edges are shown in Figure 6-2la and b. The estimation errors have 
very low sensitivity to the film element's performance. 


The performance of the three film/film-printer 

system elements were next all made off nominal at once in the 

20%, 50% and 100% amounts. The results for the twelve test edges 

are shown in Figure 6- 22a and b. The critical position estimation 

2 

error parameter, C exhibits only a 20% increase over the nominal 
value when all three spread functions at once are 100% too wide. 

A 20% increase in the variance of course corresponds to only a 
10% Increase in the standard deviation 0^ . 
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—*\ 5 M j-«~ 

20% OFF NOMINAL 
(a) 



50% OFF NOMINAL 


(b) 



100% OFF NOMINAL 
(0 


Figure 6-20 OFF NOMINAL F 1 LM/F1 LM-P R i NTE R SPREAD FUNCTIONS 
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, cr 3 2 / cr 3 2 {nominau 
FRACTIONAL EDGE LOCATION VARtANCE 


5i 


PERCENT SPREAD FUNCTION WIDTH INCREASE 
(a) 



PERCENT SPREAD FUNCTION WIDTH INCREASE 
(b) 

Figure 6-21 EFFECT OF ONE OFF NOMINAL FILM OR FILM-PRINTER 
SPREAD FUNCTION ON ESTIMATION ERRORS 


FRACTIONAL EDGE LOCATION VARIANCE . FRACTIONAL RADIANCE LEVEL VARIANCE 



PERCENT SPREAD FUNCTION WIDTH INCREASE 
(a) 



lb) 


Figure 6-22 EFFECT OF HAVING ALL FILM OR FILM-PRINTER SPREAD FUNCTIONS 
SIMULTANEOUSLY OFF NOMINAL IN EQUAL AMOUNTS 


6. 3. 4 


EBR and Film/ Film-Printer Spread Functions Off 
Nominal Simultaneously 


Tile -film and' film- printer spread functions were 
considered in combination with a 20% off nominal EBR spread function. 
The 20% value was chosen simply because it is quite gross and would 

be easy to detect but nonetheless not so gross a degradation as to be 
totally unlikely to occur. The results of the computation for 20%, 

50% and 100% width increase of a single film-printer spread function 
are shown in Figure 6-23 0 Comparison of these data with the 
preceding shows little difference. Thus to hold a reasonably coarse 
tolerance an EBR performance does not preclude the even coarser 
tolerances indicated for the film/film-printer spread functions. 


6. 3.5 


Film Development - H-D Curve 


The following target edges were generated in the 
computations discussed in this subsection. 


Designation No. 


1 

2 

3 

4 



(normalized (normalized 

radiance level) radiance level) 


(edge 

location) 


. 25 . 2 50. 

" .8 ” 

.65 .2 

" .8 » 


A 



Designation No. 


. 05 1 

. 15 4 

.45 3 

.55 2 
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OFF NOMINAL SIMULTANEOUSLY 



PERCENT SPREAD FUNCTION WIDTH INCREASE 



The processing Y of each development process was 
varied and the change in variances computed. 

First Generation 

The set of H-D curves used is shown in Figure 6-24. 
Nominally Y, = 2.0. Computations were performed for off nominal 
Y f = 1. 8, 1.9, 2.1 and 2. 2. The results are shown in Figure 6-25. 

Note: 1) The very large sensitivity to Y compared to the preceding 

spread function results. 

2) The difference in sensitivity of each target edge to the processing 
changes. The system begins to exhibit nonlinear response for small 

departures from nominal Y . 

3) Edges 2 and 4 which both have A^ = . 8 have rapidly 

increasing estimation errors for > 2. 1. The increased Y evidently 
causes saturation and clipping at this radiance level. 

4) Increasing Y actually decreases the estimation error 

to some extent before the catastrophic increase for the high radiance 
level targets occurs. 

Second Generation 

Nominal second generation processing is to = 1. 0. 

Computations were performed for = .8, .9, 1.1, and 1.2. Results 
are shown in Figure 6-26. Sens itivities are very high for all targets. 

From the cr Z sensitivities, it is clear that the high radiance values 
are being forced into saturation and are clipping. Consideration of the 
crj 1 sensitivities indicates that the only target with relative insensi- 
tivity to Y is target number one with normalized radiance values 
. 25 and . 2. Target number four which has the highest radiance level is 
most sensitive to an increase in ^ . The graphs show that it is 
critical to hold this processing element as close to the nominal value 
as possible. 
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D OUT ■ 



D IN 

Figure 6-24 FIRST GENERATION H-D CURVES 


87 


{ r, / 0^ (NOMINAL) 

FRACTIONAL RADIANCE LEVEL VAR] ANCE 



Figure 6-25 EFFECT OF OFF NOMINAL FIRST GENERATION PROCESSING 7 
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Figure 6-26 


EFFECT OF SECOND GENERATION PROCESSING 7 




Third Generation 


Nominal, third -generation processing is to ^ = 1.0. 

. O 

Computations were performed for ^ = .8, .9, 1.1 and 1.2. 
Results are shown in Figure 6-27. Sensitivities are very high for 
all targets as noted for the second generation processing. 


6. 3. 6 Measured Nonlinearity Data 

The simulation was run using the H-D curves 
measured from NDPF samples discussed in Section 6. 1 and the 
set of twelve edges defined in Ta"ble 6-2. Using the data obtained 
from the fifteen step calibration grey scale on an actual frame of 
imagery yielded the plot of edge position estimation error shown 
in Figures 6-28 and 6-29. Note that the data points lie quite close 
to the nominal results. The spread of estimation errors at each 
value indicates that the system is not linear - the estimation 
error depends on the absolute A^, A^ values as well as|A^ - A^.j 
There is no evidence that the frames examined are typical NDPF 
products. 


The twenty-one step quality control grey scale H-D 
data was also used. For this data, the system exhibited extreme 
nonlinearity which resulted in severe clipping. Figure 6-30 is a plot of edge 

position estimation error versus ^/-p . It is included to dramatize 

the effect processing control can have on the image quality. As the 21 
step grey scale data are inconsistent with both the Handbook (nominal) 
and the measured 15 step grey scale data, it is likely not typical of 
NDPF performance. 
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SOLID LINE IS CURVE OBTAINED 
USING HANDBOOK NONLINEARITY DATA 



6. 4 RBV Output Product Quality 

Since the RBV sensor has been inoperable during the 
ERTS-A mission, the RBV configuration NDPF performance was not 
as thoroughly examined as for the MSS, Referring back to Figure 6-1, 
we see that the RBV spread function is somewhat worse than the' MSS 
spread function. However, when the D/A conversion and Bessel filter 
operations are added to the MSS sensor performance, the RBV input 
signal, which is analog to begin with, has somewhat better frequency 
content than the MSS input signal. Consequently, the estimation errors' 
for the RBV products would be slightly smaller than those for the MSS 
presented above. In addition, the sensitivity of the estimation errors 
to changes in the EBR and film-printer spread function would be 
somewhat greater. But since these sensitivities are so low and the 
difference between the RBV and MSS analog inputs so small, repetition 
of the computations for the RBV was not required. It is concluded 
that quality control procedures -adequate for the MSS products are 
adequate for the RBV products as well. 

6. 5 Effect of Better Sensor Performance 

It was desired by NASA to obtain an estimate on the improve- 
ment in image quality that might be obtained with better sensor performance 
(Figure 6-31). It is emphasized that use of a better payload would require 

reexamination of two approximations valid for the current system: 

1. The assumption that remains constant with small 
changes in subsystem element performance becomes in- 
valid when the processing element spread functions are 
not narrow compared to the sensor spread function. 

2. The linearity of the contact printing process is ques- 
tionable for imagery with higher cutoff frequency than 
that imposed by current sensors. 
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Figure 6-31 HYPOTHETICAL IMPROVED SENSOR SPREAD FUNCTIONS 
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If the assumptions were to remain valid, and the sensor 
performance were improved by factors of two and four, the edge position 
estimation errors would improve as shown in Figure 6-32 . The RBV 
spread function was used to obtain these data in order to avoid having to 
change the D/A converter and Bessel filter spread functions together 

with the MSS sensor. Clearly doubling the sensor "resolution" does not 
reduce the variance by a factor of two. 

We believe that the first assumption could be invalid 
for the improved cases considered. However, the technique is still 
valid: the term would simply have to be computed as a function of 

element performance as well, which involves slightly more time, but 
involves no conceptual difficulties. If the payload spread function becomes 
close to processing element spread function sizes, it may be necessary 

t 

to use a two dimensional simulation in order to correctly, determine the 
effect of element performance on the two dimensional noise field. 

Should the contact printer linearity become invalid, a 
better mathematical model for that physical process than currently 
exists will have to be found. 

6. 6 Absolute Performance 

To relate the relative measurements considered thus far 

to absolute estimation errors requires determination of N . A reliable 

o 

estimate of the sensor noise level was not available, and, consequently, 
the following "absolute" performance calculations are for illustration only. 

It is clear from Section 3. 3 that N depends on the aperture used in the 
measurements. The long, narrow slit approximation will be used here 
for illustration. For the current ERTS system, this approximation applies 
to slits narrower than seven micrometers . 
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FRACTIONAL EDGE LOCATION VARIANCE 
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The values for Oj A ^ presented earlier in Figure 6-7 
for the "copy film" were substituted into Equation 3-21 to compute 
N q for slit lengths of 30, 150 and 300 micrometers. This corresponds 
to targets about 100, 500 and 1000 meters long on the ground. The 
resulting N q values are shown in Figure 6-33. These data were combined 
with the results for the bar target radiance level (Figure 6-15) and bar 
width (Figure 6-16) using Equation 3-2 to obtain the absolute percentage 
errors for the three lowest contrast targets. The high contrast targets 
were excluded to permit the use of an N q value at the average signal level 
(transmittance) of the target and avoid the dependence of the noise upon 
signal level. The resulting percentage errors are shown in Figure 6-34 
as a function of the width of the target in meters on the ground for a 
target 500 meters long. Both the radiance and width estimation errors 
decreases rapidly with increasing target width as would be expected. 

The independence of M ^ on target contrast, AR/Rm, combined with 
the increase of noise power with radiance or transmittance cause the 
percentage error in radiance estimation to increase rather than decrease 
with contrast or radiance level. The width estimation error, on the other 
hand, decreases with increasing contrast. The reader is reminded that 
these error estimates assume optimal data use and that their accuracy 
depends directly upon the accuracy of the system element performance 
descriptors used as input to the IDSS Program. 

6. 7 Geometric Mapping Error 

The positional estimation errors calculated 
apply to local measurements. For measurements across the frame 
format, geometric mapping error must be added. No new work was 
done in this study to assess the magnitude of this error. Residual 
errors remaining after application of correction algorithms currently 
employed by the NDPF are summarized in Appendix A. 
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WHITE NOISE POWER DENSITY -Mm 



Figure'-6-33 NOMINAL WHITE NOISE POWER DENSITY 
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% ERROR IN WIDTH ESTIMATE % ERROR IN RADIANCE ESTIMATE 



BAR TARGET WIDTH (METERS ON GROUND) 
(b) WIDTH ESTIMATION ERROR 


Figure 6-34 TYPICAL VALUES FOR PERCENTAGE ERROR IN TARGET WIDTH AND 
RADIANCE MEASUREMENT FOR A RECTANGULAR TARGET 500 
METERS IN LENGTH 
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6.8 


Radiometric Estimation Error and Calibration 


The radiance estimation errors calculated 
yield only the local error and must be added' to the large scale error. 
Algorithms for correction of non-uniform sensor sensitivity are 
currently used in the NDPF. They are summarized in Appendix A 
along with the residual errors. 

To determine the target radiance level users must 

relate output product transmission to ground radiance by removing 

the effects of the atmosphere and contamination of sensor optics. 

For the past several years we have been developing techniques to 

calibrate the transmittance of aerial photographs in terms of ground 
21-23 

reflectance. We examined the application of these techniques 

to some typical ERTS imagery and found that they might be successfully 
applied. 


The calibration techniques use analyses of scene shadow 
areas to establish the relationship between sensor irradlance and 
ground reflectance. Very briefly, a linear relationship exists between 
the radiance just inside and outside a shadow. The parameters of this 
relationship are related to the component of radiance due to atmospheric 
scattering or flare and to the ratio of sunlight to skylight irradiance. 
Thus, the relationship between sensor irradiance and scene radiance 
or reflectance can be established by measuring sunlight and shadow 
radiance for a set of shadows and determining the parameters of the 
fit to the data set. The shadow calibration techniques require no 
atmospheric modeling (and consequent model parameter measurement); 
rather, the calibration proceeds entirely from the sensor record of 
shadow elements within the scene to determine the atmospheric 
component of radiance. 

As a result, a calibrated sensor record is obtained 
in which sensor response is directly proportional to scene reflectance. 
Absolute calibration requires knowledge of only one element in the scene, 

U.S. Patent 3 , 849, 006. 
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similar to laboratory use of a MgO standard or its equivalent. It 
should be noted that the calibration process can be automated so that 
data tapes are immediately corrected for the effects. 

. . Our previous analyses have been at sufficiently large 

image scales that shadows from buildings and boulders could be 
utilized. The scale and resolution of ERTS imagery preclude use of 
such shadows; however, we have found that shadows on terrain with 
large relief and certain cloud shadows can be used to calibrate ERTS 
imagery successfully. 

The first step in the calibration procedure is to reduce 
the measured film densities or sensor voltages to relative exposure 
or radiance by means of a D-log E or instrument response curve. 

Analyses then proceed from these relative exposure values. It is 
important to note that it is not necessary to have a relationship between 
sensor response and absolute exposure. All analyses can proceed 
from a relative exposure relationship; hence, the problems currently 
existing with the MSS mirrors (contamination) have no effect on the 
calibration process. 

The resultant exposures are dependent upon meteorological 
conditions, altitude of measurement and illumination conditions such 
as proportion of sunlight to skylight and the amount of "airtight" (the 
•contribution to exposure by illumination scattered to the sensor by the 
atmosphere beneath the satellite). These effects are depicted in Figure 6-35. 

All of these effects can be approximately coupled into 
three parameters for a given spectral band: <x , a.' t and . The 

parameter <X is proportional to atmospheric transmittance and total 
(sunlight + skylight) irradiance; cc is proportional to atmospheric 
transmittance and skylight irradiance; and is proportional to the 
amount of air light in the scene. The exposure in sunlight, E, for 
an object with reflectance R is 
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Figure 6-35 ATMOSPHERIC AND ILLUMINATION EFFECTS INVOLVED IN 
ESTABLISHING EXPOSURE OF A TERRAIN ELEMENT 


10 


E - <XR + ft , 


(6-1) 


while the exposure of the same object in shadow, E ! , is 

E' = cL R + ft (6-2) 

Measurements of terrain reflectance thus requires 
knowledge of a , tx and ft in each sensor spectral band. These 
parameters can be determined in a straightforward manner using a 
shadow calibration procedure called the Scene color Standard (SCS) 
technique. 


Calibration is accomplished by densitometry of the 

21 

illumination discontinuities at shadow edges. It is convenient to 

/ 

write ix - 6 + <x with 6 a term proportional to solar irradiance 

only. The discontinuity measured at shadow edges on two different 
terrain elements then determines ft and 6 / cl as follows. 


In the sunlight just outside a shadow the exposure, E, 
is 

E = ( cf + cxl )R + ft (6-3) 

The R is the terrain reflectance. Just inside the shadow the exposure 
E r is given by Eq. (6-2), Eqs. (6-2) and (6-3) yield 


E = (1 + — ) E 1 - ft A. (6-4) 

a. ck 

Eq. (6-4) is a linear relationship between E and E' 
d cf 

with slope (1 + —pr ) and intercept -ft . Two similarly oriented 

shadows determine the slope and intercept, and hence ft and -p- t 
In practice a number of shadows are analyzed, and a least squares fit 
is made to the data. 


The atmospheric conditions have now been determined. 
One aspect remains -- that of establishing an absolute level of 
reflectance, akin to laboratory use of a MgO standard or its equivalent. 
A tar or sheet asphalt scene element in sunlight (roadway, roof) is 
usually used to establish the value of cc and complete the calibration. 


ORIGINAL PAGE IS 
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These elements are used as: (1) their reflectances are spectrally flat; 

(2) their reflectance remains constant .to good approximation .ove-r -the - - 
- - - ’year;' and "(3j”their reflectance can be easily, estimated or measured. 

Other objects more appropriate for a particular image can, of course, 

.be used.. 

The SCS technique was applied to several- frames of 
ERTS-A imagery. An example of calibration using terrain shadows 
is given' in Part a of .Figure 6-36, a print .of Band 6 (IR) ERTS MSS 
image of Antarctica taken 14 February 1973. Some typical shadows 
used for calibration are depicted by arrows on the figure. About 
twelve shadows were used in ihe calibration, although easily fifty 
shadows could be found and so utilized.. The reflectance of the snow 
was established at 70% and used as the standard in this case. 

Figure 6-37 displays the linear fit to sunlight and 
shadow exposures for the Antarctica scene. The parameters beta and 
illumination ratio show that the atmospheric component of radiance 
at the satellite is equivalent to a 7% scene - reflector-, and that the ratio 
of sunlight to skylight irradiance is 3.2. Equivalent atmospheric flare 
of 7% means that a perfectly black object (zero reflectance) would appear 
to be a 7% reflector, a true 5% reflector would appear as a 12% reflector, 
etc. 


The snow covered areas in the scene are regions of 
70% reflectance. Even for these bright areas, 10% of the scene radiance is 
due to atmospheric flare. The rocky peaks are measured from the 
ERTS image to have a 17% reflectance so that over 40% of the apparent 
radiance is due to atmospheric flare. ERTS investigators using the 
multispectral characteristics of scene objects clearly must take great 
care to remove the effects of atmospheric flare from their data before 
analyses. 
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(a) BAND 6 - ANTARTICA 
(TERRAIN SHADOWS) 


(b) BAND 6 - LAKE ONTARIO 
(CLOUD SHADOWS) 


Figure 6-36 


ERTS MSS IMAGERY USED TO DEMONSTRATE RADIOMETRIC CALIBRATION 


Figure 6-37 



RELATIONSHIP BETWEEN SENSOR RADIANCE IN SUN AND SHADOW 
FOR BAND 6 OF ANTARTICA SCENE 
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RELATIVE EXPOSURE IN SUNLIGHT 



RELATIVE EXPOSURE IN SHADOW 


Figure 6-33 RELATIONSHIP BETWEEN SENSOR RADIANCE IN SUN AND SHADOW FOR 
BAND 4 OF LAKE ONTARIO SCENE 
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SECTION 7 


CONCLUSIONS AND RECOMMENDATIONS 

.In order to predict image data system performance 
and evaluate product quality a new technique was developed, implemented 
and demonstrated using the NASA Data Processing Facility (NDPF). 

The technique has wide potential application, not only in quality control, 
but in system design and performance prediction. 

Quality was defined as the limiting accuracy in estimating 
information from an output product in terms of the Cramer-Rao error 
bound. An Image Data System Simulation (IDSS) software program was 
written to compute this error bound in terms of accepted system element 
performance descriptors. The development and implementation of this 
technique are major results of the study. Insights into quality, control 
requirements for the NDPF are also among the significant results. 

7. 1 NDPF Quality Control Requirements 

As a result of the application of the IDSS program to 
the NDPF requirements for quality control were reviewed. These 
requirements are discussed below in terms of OTF monitoring and 
photographic processing control 

■7. 1. 1 OTF Monitoring 

The sensitivities of estimation errors to changes in 
spread function widths have been shown to be quite small. Changes 
of fifty percent are possibly tolerable, changes of twenty percent are 
certainly tolerable. Since the increase in spread function width is 
manifested in frequency space as an equal fractional decrease in OTF 
cutoff frequency, rough methods for detecting changes in cutoff fre- 
quency are adequate. Resolution targets, accepting the variability 
inherent in the human element, can be used and are recommended 
because of the speed of diagnosis and relatively simple operator task 
involved. 
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In addition, an edge gradient spectrum computation pro- 
gram has been written (Appendix D) and can be used in diagnosis. 

Edge target signals can-be- generated electronically and used as input 
to the EBR. These targets can be measured on an output (3rd' generation) 
product to define .the quasi-linear region of signal input as demonstrated 
in Section 6.2. Under nominal conditions the NDPF is linear over a 
signal range at least from 0.2 tp 0.8 normalized radiance. Measurements 
of these targets on intermediate products will permit the assessment of 
the frequency degradation for each stage. 

If the system is constrained to nominal (Handbook) per- 
formance, the edge gradient procedure can also be applied to traces of 
ground- objects to yield the overall system frequency response. When 
the system behaves linearly, changes in the sensor OTF can thus be 
remotely monitored. 

Most of the spatial frequency or resolution degradation 
occurs prior to production of the archival copy, i. e. at the spacecraft 
sensors. Improvement in sensor OTF will produce an improvement in 
product quality. With improved sensor performance, the need to monitor 
system OTF will become more significant. 

7. 1 . 2 Photographic Processing Control 

The simulation results showed the estimation errors to 
be highly sensitive to changes in H-D curve If . Deviations of less than 
10% in 'if decreased product "quality" by a factor of two. The shape 
of the individual processing element H-D curves must therefore be as 
tightly controlled as is possible. Current quality control procedures 
measure the H-D curves and routinely perform a polynominal fit to the 
measured data. However, the entire curve is not observed, rather the 
values of only certain selected- points are checked. 
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It is recommended that the complete carve be monitored. 

A graphics output could be attached to the computer which' cur rently 
performs the polynomial fit, or if that is not feasible, a programmable 

desk calculator with graphic output could be substituted for the computer 
irx this application. In any event,- the r-esults show that it is essential 
that the entire H-D curve be constrained to remain within tight bounds 
if the loss in quality and departure from linearity is to be avoided. 

It would also be useful to monitor composite ; that is, 
to plot third generation density against the log of the EBR input signal. 
Measurements could be made on the fifteen step calibration grey scale. 
The composite If should, of course, equal minus one, and the shape of 
the curve should be very nearly a straight line over the input signal range. 
Since the input signal that generates the grey scale does not change with 
time, it would only be necessary to measure third generation density. 
Plotting and mathematical processing could easily be accomplished with 
a programmable desk calculator equipped with a plotter attachment. 

7. 2 Potential Radiometric Calibration Procedure 

Output product transmission cannot be related to ground 
reflectance correctly unless atmospheric and payload optics contamina- 
tion scattering effects can be compensated. The SCS technique described 
in Section 6. 8 could be applied, either by the user, or perhaps at the 
NDPF to obtain relative reflectance values from transmission measure- 
ments. A reflectance standard on the ground would be required to relate 
the relative scale thus obtained to absolute values. The precision with 
which this technique can be applied to ERTS imagery has not been fully 
established, and additional empirical work would be required before an 
estimate of the error in ground reflectance can be given. 
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Comments on User "Optimality" 


7. 3 

The variances given by the jCr-ame-r—R-ao bound as applied 
to the ERTS NDPF quality control problem characterize the estimation 
error that would be obtained by using an optimal information extraction 
procedure. It is likely that the errors actually obtained by most users 
will not be as small as these. 

The term "optimal" has been used in a strict communica- 
tion theoretical sense. That the user's equipment is working properly, 
he uses the calibration grey scale properly, etc. is taken for granted. 

It is implicit that the user has a priori knowledge of the type of object 
he intends to observe and that he has knowledge of the system performance, 
and noise spectrum. Unless his detection technique makes use of this 
knowledge, the estimation errors predicted by the Cramer-Rao bound 
will not be achieved. 

7. 4 NDPE Nominal Product Quality 

The IDSS program was used- to evaluate the nominal 
quality of the NDPF 3rd generation imagery, a primary ERTS product. 

The quality assessment summarized here depends on the absolute accuracy 
of the data base supplied as hypothesis by NASA. Although we cannot 
stipulate the accuracy of that data, we state the product quality deduced for 
completeness. Quality was measured in terms of the limiting errors in 
estimating target apparent radiance, size and location. The percentage 
error in both the metric and radiometric measurements exhibited the 
expected variation with target size showing a rapid increase for widths 
less than 120 meters. For widths greater than 120 meters the target size 
measurement errors varied from 2% to 30% and target apparent radiance 
measurement errors from 0. 5% to 3. 5% depending upon target size and 
contrast. The radiometric accuracy due to large scale effects (e.g. EBR 
non-uniformity, printer shading, film sensitivity variations, ecc. ) is 
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stated to be 5% in the Data Users Handbook. Thus the ability to measure the 
transmittance of the photographic product and knowledge of the transmittance 
to. apparent radiance calibration curve have about .the same magnitude errors. 
For the larger targets the imprecise knowledge of the calibration curve 
could limit the radiometric accuracy if the user employs near "optimum" 
measurement techniques. 


7. 5 Recommendations 

Asa result of this study the following recommendations 

are made. 


(1) More complete monitoring of the NDPF film 
processing functions should be developed. 

(2) The use of the shadow- sunlight measurement 
technique for radiometrically calibrating ERTS imagery should be 
studied to establish its utility and accuracy. 

{3} Application of the technique developed under this 
study to other image data system design/analysis problems ought to 

r 

be considered. 
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APPENDIX A 

BACKGROUND INFORMATION 


To provide background information for this study a 
number of surveys were performed at its onset. These included a 
detailed survey of the NASA Data Processing Facility (NDPF) operations 
to identify the image data flow, radiometric and geometric correction 
procedures, element performance specifications and quality control 
procedures. A survey of the state-of-the-art in image quality analysis 
and evaluation was made to identify potential measures of product 
quality and element performance descriptors. Finally a review 
of users' requirements for product quality was made. The information 
obtained form these surveys is contained in this appendix and summarized 
in the main body of this report. 

A. 1 NDPF DESCRIPTION 

The information on NDPF operations was obtained from 
on-site visits, discussions with NASA personnel, NASA furnished 
documentation and the Data Users Handbook, 


A. 1. 1 


SUBSYSTEM ELEMENTS AND DATA FLOW 

The initial task undertaken was identification of the 


NDPF subsystem elements and data flow. The only information source 

available initially was the NASA ERTS Data User's Handbook . Sub- 

2 

sequently, a conference was held at Goddard SFC for the purpose of 

obtaining fundamental information for the study. In addition to the verbal 

3 -24 

information, copies of a number of NASA documents and Bendix Corp. 

25-34 

technical memos were obtained. The information contained in some 

of these documents was not necessarily current. The User's Handbook 1 
on the other hand is kept up to date; consequently, that documents plus 
verbal information received at the technical conference served as the major 
quantitative information source during the first phase of this study. However, 
some of the data contained in the Handbook is clearly idealized. The data 
were ultimately supplemented by measurements from ERTS data products . 

Note: Superscript numbers refer to the References in Section A . 4. 
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This survey was explicitly limited to consideration of the 
NDPF video to photographic output product conversion process. NDPF 
digital output products were excluded. The input to the NDPF is video 
tape recordings of telemetry from the ERTS payload. Payload and 
telemetry link elements, are external to the NDPF. However, characterization 
of the video tape input will be found to enter the analysis of the NDPF per se 
and is given a priori. 

A block diagram of the NDPF data flow is given in 
Figures A-l to A-3. This block diagram served as a basis for system 
modeling. The following comments refer to the circled numbers in 
Figures A-l to A-3. 

The Return Beam Vidicon (RBV) ( 7 ) consists of three 
boresighted RCA vidicon systems in three different spectral bands. 

Radiometric and geometric calibration capability exists in the payload. 

Radiometric Calibration Images (RCI's) consisting of white, 
gray, and black constant field illumination are taken before and after 
each sequence of exposures (10-15 minute period). These images are 
processed in the NDPF as if they were real scenes and consequently yield 
a radiance map that can be fed back into the system for updating the 
radiometric correction. Currently, the NDPF assumes that the camera 
radiometric response is slowly varying, hence the feedback or updating 
is applied to the processing of subsequent ground image sets, not to 
processing of the set corresponding in time to the' particular calibration 
images,. 


Geometric correction is provided by analysis of a reseau 
fixed at each of the vidicon faceplates. The reseaus allow correction 
for geometric optical distortion and boresight errors. A record is 
being kept of individual camera relative position and distortion. 
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Figure A-2 NDPF DATA FLOW BLOCK DIAGRAM (PART 2) 
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The RBY telemetry (z) is analog with 3. 5 mHz bandwidth. 
The signal to noise ratio on the videotape record is 12-13 db (Ref. 2). 

The three channels are expo.sed. simultaneously, but broadcast sequentially. 
The band limit of the RBV videotape record is given in Figure A-4 (Ref. 1, 
p. A-5). 


The Hughes Multi spectraL Scanner (MSS) (5) has four 

conjugate linear detector arrays, each in one spectral band. No 

geometric distortion is introduced by the instrument between the 

individual bands (registration), however a mapping error can be 

introduced. Two procedures take place in the payload for the purpose 

of radiometric calibration; a lamp provides a linear intensity wedge, 

and a pinhole image of the sun is passed over the array. The MSS 

telemetry (4) is digital with 16 MHz bandwidth. The error rate 

5 

on the videotape is 1/10 (Ref. 2). The band limit of the MSS video- 
tape record is given in Figure A-5 (Ref. 1 p. A- 13). MSS data is 
D/A converted ( 5 ) in the playback operation^ Radiometric corrections 
are applied in the digital domain before conversion. 

The Bulk Image Annotation Tape (BIAT) (6) generates 
tick marks, denotes the scene, and provides input for geometric correction 
for payload attitude. Attitude data are provided for nine points north to 
south across each frame. The frame format is divided by a nine by nine 
grid, and attitude correction is linear over each sector (only translation 

and rotation operations line by line within each sector). Attitude correction 
is the only geometric correction applied to MSS data. Attitude correction 
for RBY data is only part of the total geometric correction. 
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Figure A-5 MEASURED MSS SYSTEM SQUARE WAVE RESPONSE 




The Electron Beam Recorder Image Correction (EBRIC) 
tape ( 7 ) , provides input for radiometric and geometric correction of 
RBV imagery. 


Radiometric correction data are derived from the Radiometric 
Calibration Images (RCIs) described in the comments on the RBV, The 
Bulk processed RCIs are taken to the Precision Processing Subsystem 
where EBRIC radiometric correction tapes are generated. These EBRIC 
tapes are subsequently applied to new Bulk imagery. 


Geometric correction results from checking the 
reproduced reseaus in Precision. Comparison of the actual reseau 
coordinates with a priori knowledge results in generation of geometric 
correction coefficients for the EBRIC tape. 

The only effect EBRIC tapes have on MSS imagery is 
correction for errors internal to the NDPF. 

The Electron Beam Recorder (EBR) (If) produces all 
archival latent images. The film type used is Kodak SO-438 (Ref. 2). 

A fifteen step gray scale is put on each image. 

The first processor ( 9 ) is a Kodak Versamat used only 
for processing archival images. 

The quality control blocks (QC) (fo) consist of standard 
Kodak chemistry quality control plus the placing of a special target on the 
head and tail of each roll processed, Tbe target consists of a gray scale, 
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five equal constant density patches over the format, and a standard Air 
Force bar target in the frame center. The gray scale is read and a 
Hurter-Driffield (H-D) curve fit to. the data points. The density values 
at two exposure levels are plotted and deviation from nominal values 
used as a processing quality criterion. The constant density is read at 
the five format positions to provide uniformity data. The Air Force 
tri-bar target allows determination of on axis resolution. 

The Bulk Enlarger (fl) is a modified Miller-Halzwarth 

The Strip Printers (12) are Kodak Colorado printers. 

The processors (L3) are Kodak Versamats and the paper 
processor (14) is likewise a standard Kodak device. 

The selective printing process (15) is a hand 
operation. Three Mark III printers are employed. 


Registration (l£) of color composites is accomplished 
by punching holes in the black and white transparencies and fitting these 
holes on locating pegs in the spectial Bendix Color Composite Printer (l?) 
Registration accuracy is one pixel. 

The color processing is done in a Kodak 1811 Versamat (18) 

Ground control points (l?) are provided by film chips cut 
from a contact print of a master plate. 
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The Viewer /Scanner, Video Processor, and Scanner 
Printer (20) are analog devices. The Viewer/Scanner image is 
digitized and the correlation (21) and fit takes place in the digital 
domain. The Precision Processing Subsystem generates either 
precision images or EBRIC tapes. No radiometric correction is 
performed on products in precision. The system is expected to meet 
geometric specifications but not radiometric specifications (due to 
the optical transfer functions of the various devices). Accuracy 
within each precision processed frame is two hundred feet with 
respect to ground control points. Absolute ground location accuracy 
is not specified. 

The annotation tapes (24) provide tick marks and annotation 
for Universal Transverse Mercator (UTM) projection. 


Input to precision results from cutting out selected 
•second generation negatives, pasting them to a plastic sheet, and contact 
printing to produce a third generation positive. 

The photographic processing subsystems have a . Q3AD 
capability (AD = photographic density error) but this precision inhibits 

throughput. Processing to 2f'=-1 is sought for the system as a whole 
(including the EBR as well). All film (except in the EBR) is Kodak 
SO-467. Chemistry is Kodak 641 or 2420. 

A. 1. 2 BULK PRODUCT GEOMETRIC AND RADIOMETRIC 

CORREC TIONS 

The ERTS Data Users Handbook (Ref. 1) discusses 
system performance including output product geometric and radiometric 
accuracy. Two types of geometric accuracy are quantified; mapping 
accuracy and spatial registration accuracy. The error sources are 
classified into three categories; errors external to the sensor, internal 
sensor errors and processing errors. Table A-l summarizes the 
external and internal errors for both the RBV and MSS sensors. These 
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TABLE A-l GEOMETRIC ERROR SOURCES 


- 

rms' magnitude (meters)] 

ERROR 

RBV 

MSS 

External 

723 

710 

Internal 

1355 

26 

Net-NDPF Input 

1536 

710 

Bulk Product-Mapping 

766 

743 

Bulk Product-Registration 

336 

159 


combine to yield the net error in the inpat signal to the NDPF. During 
' the input videq-to-film conversion corrections are applied to reduce 
the resulting error in the bulk image products. These corrections 
are partially negated by degradations introduced by the NDPF. In 
addition, the NDPF does not attempt to correct for all of the input error 
sources; specifically external sources including exposure time 
spacecraft ephemeris and altitude errors are not corrected. Figure A-6 
presents schematics of the image correction procedures for both the 
RBV and MSS sensors. The net mapping and registration rms errors 
in the bulk film products after correction are also included in Table A-l. 

RBV radiometric calibration is accomplished through the 
EBRIC tape which contains the EBR beam adjustment information. These 
data are derived from two sources: 1) the pre-flight RBV radiance mapping 
and 2) the on-board erase lamps which produce in-flight Radiometric 
Calibration Images (RCI's). The radiometric error in bulk processed 
RBV imagery is 9% of full scale sensor voltage. The equivalent error 
in radiance varies with level because RBV voltage is not linearly 
proportional to apparent radiance. 
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Figure A-6 BULK PROCESSING IMAGE CORRECTION SCHEMATIC 





















MSS radiometric calibration is accomplished through an 
internal radiance source and sun scans. The in-flight internal calibration 
data are obtained every MSS scan and -it is us'e'd to obtain the optimum 
linear relationship between radiance in and voltage out for each detector. 

The internal calibration source was initially measured on the ground. 

Long term drift in the absolute radiance calibration is corrected through 
periodic sun scans. The MSS radiometric error in bulk processed MSS 
imagery is 5% of full scale sensor count or apparent radiance. 

A. 2 _SURVEY O F GEN ERAL USERS' REQUIREMENTS 

That the image quality characterization resulting from this 
study contain sufficient information for users to determine the adequacy 
of ERTS imagery for individual needs is a necessary condition. Users' 
ability to define "image quality" is not established and completeness 
of an image quality characterization based on users requirements would 
hot be expected. But a general understanding, at' least, of what tasks 
users of ERTS Imagery would like to accomplish is certainly required if the 
"quality" which is controlled is to have relevance. 

A literature search was consequently conducted. A 
large variety of journal articles and symposium proceedings were sampled 
to survey aerial and satellite imagery users’ activities and needs. (Those 
not referenced in the text are included in the bibliography. ) In many 
instances the user’s requirements were stated or implied only in 
general terms. Other cases provided specific, but qualitative, definitions 
of image' quality parameters. Finally, a few authors stated specific, 
quantitative, image quality requirements. 

I 

It should be noted that some of the papers concern 
aerial remote sensing, with its inherent capability for image resolution 
superior to that of ERTS imagery; thus they cannot be directly applied to 
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ERTS data user requirements. They do, however, provide additional 
background on the use of remote sensing in the field discussed. 


The surveyed articles were grouped into broad classes: 
land use studies, vegetation and crop studies, geology and geography, 
oceanography and hydrology, and oil pollution surveys. User requirements 
as deduced from the articles listed in the bibliography will be discussed class 
by class. 


A. 2. 1 Land Use Studies 


More articles on land use surveys were available than 
on any of the other topics. In addition, this class of articles contained 
some of the most specific and quantitative statements of user image 

35 

quality requirements. In one article concerning a study performed to 
determine the effects of various environments on the levels and types 
of information retrieved from orbital-acquired (Gemini and Apollo) 
imagery, the authors employed a technique of placing- an artificial grid 
of resolution cells over the image and then counting the number of image 
elements in that cell. They concluded that for land use mapping the 
minimum resolution requirement is that at least 50% of the cells contain 
only one element. In terms of ground resolved distance (GRD) they 
stated: 


Level of interpretation GRD 

General identification of terrain 300 ft. 

Precise identification 15 ft. 

Description 5 ft. 


Another study with a similar purpose discussed space photographs 
simulated by artifically reducing the resolution of aerial photographs 
by various degrees. The photos were then subjected to human photo- 
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interpretation. The results indicated that simulated ERTS data (GRD of 

250 ft.) contain sufficient information to allow an interpreter to discriminate 

among woody vegetation,, g.r-as-sla-nd and wafer bodies. For more detailed 

information such as species identification of the woody vegetation, the 

imagery must have a ground resolved distance of 50 feet'which is beyond 

37 

current ERTS capability. A third study supports the above results 
with its own conclusion that ERTS is most useful for broad land use mapping 
in regions such as the Great Plains v/hich is dominated by large blocks 
of natural categories. Additional articles in this area discuss high 
resolution aerial imagery or land use classification schemes. 


A. 2. 2 Geology and Geography 


Another area for potential ERTS data users is that of 

geology and geography. An example of a remote sensing project where 

both the radiance -and size of a small object is required in given by 

3 8 

Vincent and Thomson . Silicate concentration in rocks is estimated 

to within 14% using a technique that relies on the SiO^ concentration 

dependence of the shape of the emissivity spectrum of silicate rocks. 

Radiometric fidelity of lj J. in wavelength is required. A resolution 

cell of less than 100 ft. is preferred, but the authors feel this 

requirement can be reduced through the use of visible and near IR 

39 

sensors. Another paper compares aerial and space- acquired (Gemini IV) 
imagery for use in geological mapping. The authors concluded that the 
space photography identification performance was good, but not equal to 
that of aerial photography because of reduced resolution on the space 
photographs. 
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A. 2. 3 Oceanography and Hydrology 

40 

Stevenson gives a list of specific qualitative image 

derived parameters as well as some quantitative estimates of the desired 

precision in these parameters. Cited as of interest to oceanographers 

are: ocean color, sea-surface roughness, sea-surface temperature, slope 

of the ocean surface and of significant waves, atmospheric profiles of 

temperature, moisture and CO and lunar magnitude of tide producing 

^ o 

forces. Estimates of accuracy required in image quality are 100A in 

o 2 

wavlength; 1 C in temperatures; spatial resolutions of 10km for islands, 

.2 2 

coasts, and current boundaries, 500km for open oceans, 25km for ocean 
surface patterns, 5m for wave heights, and 100m for surface features 
and required repetition intervals of 24 hrs. near coasts, 5 days in the 
open ocean. 


Hydrological data goals as given by Eock were stated 
in broad terms as: study of the hydrological cycle; ice and snow, saline, 

water pollution mapping; and coastal surveys. 


A . 2. 4 Oil Pollution Surveys 

Another field for potential users of ERTS data is oil 

pollution monitoring. Several articles have implications for the 

42 

relevance of ERTS data to this field. Wobber suggest satellite coverage 
for use in monitoring open sea slicks, major coastal spills, and highly 
probable spill areas (pipelines, ports). Aerial coverage, however, is 
preferred for its speed and flexibility in covering a sudden oil spill. 

In addition, this and other articles place the useful region of the EM 
spectrum for oil pollution monitoring at visible blue to ultra-violet. 
Thermal infrared can also be useful if ground truth data are available. 
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A, 2. 5 Crops and Vegetation Monitoring 


Crop disease surveys and vegetation mapping .a r-e anothe-r- 

a-r-ea where use of aerial remote sensing makes it a candidate for an ERTS 

43 

data user, Philpotts discuss crop disease patterns revealed by moderate 
scale color IR photos, 

A , 2. 6 Summary 

Based on the literature sampled, no definition, of '‘image quality" 
useful in accomplishing the study objectives can be drawn from the users. 
Surely, one might adopt the terms "radiometric fidelity", "geometric 
fidelity", "resolution", but these terms have different meanings to 
different people and are certainly not sufficiently well defined to provide 
a useful basis for quality control criteria. They are general terms 
which classify rather than specify the ability to make certain measurements 
on photographs. 

One can look at the users' task s just described once 
again however, and ask exactly what sort of measurements is each 
user making. The following seem to span the tasks: 

1. How accurately can a boundary between different 
transmission levels be located on a photograph? 

2. How well can the radiance, siz$. and location of 
small objects be measured? 

3. How well can the separation between two objects 
or boundaries on a photograph be measured? 

4. Kow well does that separation represent the 
separation on the earth? 

5. How well can the transmission of a photograph be 

measured? 
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6. How is that transmission related to radiance at 


the earth? 


Clearly, any working definition of "image quality" adequate 
for the task at hand must be capable of obtaining quantitative answers to 
such questions and must relate those answers to measurable properties 
of elements of the NDPF, 

A. 3 SYSTEM ELEMENT PERFORMANCE CHARACTERIZATION AND . 

"IMAGE QUALITY" 

This section is to delineate techniques commonly employed 
to characterize image processing system element performance and to 
predict the quality of the images produced. The purpose is to establish 
some concepts that will subsequently be used and to point out why some 
others are unsatisfactory for the study objective. 

In order to ensure familiarity with the current state of 
the art, a literature search covering the period 1968 - date was under- 
taken. This section reflects the content of the resulting bibliography. 

Four overlapping categories, optical transfer function, 
resolution, noise, and subjective image quality, provide convenient 
headings for discussion. 

A. 3, 1 Optical Transfer Function 

The optical transfer function has been shown to be a 
useful tool to characterize the performance of many image processing 
devices. Some techniques for its measurement are pointed out in this 
section. 
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The blur introduced by an optical system linear in 
intensity can be characterized by a ''point spread function" S (X) defined 
implicitly by the convolution:. 


o ( iT) s c 5) 

where: i(X) = image brightness, OCX) = object brightness, 
a position coordinates in plane orthogonal to optical axis, 
convolution). 


(A-l) 

and 

(■&■ denotes 


Thus: . (A -2) 

o (T; t- ( 

where: X (V)^£ (~x) , O(t>)x=*0(x) , T(l?)l=Z S(X) , and 

denotes a Fourier transform. 


the optical transfer function (OTF). It is in general a two- 
dimensional, complex valued function, j T (V) | is the modulation » 
transfer function (MTF). It is emphasized that the OTF is not a 
measure of .image quality but merely the frequency response function 
of a linear device and consequently a measurable property of the per- 
formance of that device. 

Application of the OTF concept to photo-optical systems 
requires linearization of the generally non-linear development "process. 
Generally the OTF concept is applied to the object to exposure image 
transfer process. 


Methods of measurement of the OTF of photo -optical 

44-47 

systems are well-documented in the literature. Such measure- 

ment obviously requires use of object targets with known Fourier spectra. 
The three most commonly employed are considered here in brief. 
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1) Sine Wave Targets - Sine wave targets consist of 
objects whose brightness varies sinusoidally with distance (Figure A -7). 
The targets are photographed, image transmission is measured and 
transmission values are converted to exposure space through the Hurter- 
Driffield (H-D) curve. The ratio of the exposure contrast to the target 
contrast is the value of the MTF at the target spatial frequency. The 
phase component of the OTF cannot be measured unless a method of 
measuring the displacement of the image sine wave peaks from the optic 
axis (or other convenient datum) is provided. 


TARGET 



tM 


IMAGE CONTRAST 
TARGET CONTRAST 



2) Square Wave Targets - Square wave targets consist 

of objects whose brightness consists of a periodic discontinuous variation 

between brightness levels (Figure A -8a) . When the widths of the two 

levels are the same, the targets can be used in the same manner as 

sine waves. The "square wave modulation" obtained is not the OTF. 

A relationship between square wave modulation and the modulus of the 
48 

OTF is: 
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Figure A-8 SQUARE WAVE TARGETS 



ro>) 

where M 0 is the target contrast and 3^9) is the image square wave 
contrast. 

Another use of square wave targets results from observing 
that a target consisting of a finite number of square bars of arbitrary- 
width and spacing can be represented by an infinite periodic array of 
Dirac S functions convolved with a square pulse equal to the bar 
width and multiplied by a square pulse equal to the target width (Figure 
A-8b). The target spectrum is consequently given by an infinite periodic 
array of & functions convolved with the Sine ) function corres- 

ponding to the target width and multiplied by the Sine function corres- 
ponding to the bar width. The photographic image can be scanned 
with a micr odensitometer , converted through the H-D curve to exposure 
space and Fourier transformed. The ratio of the calculated exposure 
image spectrum to the known target spectrum is the OTF, 

3) Edge Targets - Knife edge targets consist of a step 
function in object brightness (Figure A -9). Since the derivative of a 
step function is a <5 function, the spectrum of the target gradient is a 
constant. Consequently, the OTF is given by the Fourier transform 
of the derivative of the exposure image of an edge obtained by 
scanning the photograph with a microdensitometer and using the 
H-D curve for conversion to exposure space. Edge targets have found 
wide acceptance because they are easy to produce, analysis, is straight- 
forward, and they have the additional attractive property of often being 
found in natural aerial photographic scenery. 

It should be noted that the targets discussed all yield 
one-dimensional transfer functions, 

ggMSS 
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Figure A-9 EDGE TARGETS 

The accuracy of OTF measurement in photo-optical 
systems is related to system photographic (grain) noise. Four inde- 
pendent studies of the relationship of noise to OTF accuracy are refer- 
44 49-51 

enced. ’ All concern delination of proper filtering techniques 

to minimize the uncertainty introduced by grain noise. The work of 

44 

Kinzly and Mazurowski deserves special mention since, in addition 

49 

to application of the technique suggested by Blackman , it develops 
a promising adaptive filtering technique. All four references have 
primary application to edge gradient techniques. In consideration of 
optimum methods for reducing error due to noisy data, the obvious 
should not be overlooked: since the measurement techniques are one- 

dimensional, noise in the data can be reduced by increasing the length 
of target and scanning (slit) aperture in the direction orthogonal to the 
scan. 
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Two image quality assessment parameters based on 

the OTF are the Strehl definition and Shade's equivalent passband. 

The Strehl definition is the ratio of the volume enclosed by the measured 

MTF to the volume enclosed by the diffraction limited MTF. Shade's 

, .2 

equivalent passband is the volume enclosed by | , Roet-ling, 

52 

et al. have shown that Shade’s equivalent passband is equivalent to 
acutance (mean square edge density gradient) and is consequently use- 
ful as an expression of detail rendition as perceived by a human observer. 

A. 3. 2 Resolution 

Resolution is the ability to distinguish between two adjacent 
objects in an image. A number of resolution criteria are in use. 


Rayleigh's criterion is that the diffraction limited images 
of two points are just resolved if the central maximum of one lies on 
the first minimum of the Airy disc of the other. Thus, Rayleigh's 
criterion can be expressed: 


£ = 


0.61 A 

N A 


where: X = wavelength 

NA = numerical aperture 


(A - 4) 


A perceptible dip exists between point images separated by Rayleigh's 
criterion. The dip disappears at 

*- °-wr ' (A - 5) 

which is Sparrow's criterion. 

Both the Rayleigh and Sparrow criteria are clearly 
related to the spread function and can thus be derived if the OTF is 
known. 
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A number of resolution criteria are simply defined by 
an observer's ability to distinguish the existence of a particular target. 
Such criteria depend not onl-y -on the properties of the image but on the 
properties of the detection process as well. The most common is the 
standard Air Force Tri-Bar target. One can obtain a "modulation 
detectability curve" by having a number of subjects observe tribars 
of varying spacing and contrast and plotting the statistical detection 
threshold contrasts versus spatial frequency, "Resolution" can then 
be defined as the intersection of an MTF and a modulation detectability 
curve. Uncontrolled variables and experience produce uncertainty 

in this measure of resolution. One summary measure of "image 
quality" that is in use is the area enclosed between the modulation 
detectability curve and the MTF, 

The term "resolution" is sometimes applied to the ability 
of an optical system to "resolve" a specified object. This definition is 
similar to the preceding one but requires recognition as well as detection. 


The concept "resolution" can be more precisely 
expressed as estimation error as a function of object properties 
and imaging system performance. In this case we eliminate the 
subjective detection process in order to obtain increased predictive 
validity. 

A . 3 , 3 Noise 

In photo-optical systems, the major noise source is the 
granularity of the emulsion and is expressed in the granularity constant, 
G. The rms density fluctuation observed in scanning a uniform density 
area is 



where: A = area of scanning aperture 


(A-6) 


154 



In actuality the emulsion records the continuous exposure distribution 
as a discrete, thin but nevertheless three-dimensional, distribution 
of silver particles. The photographic macro-image is a continuous 
intensity distribution which results from multiple scattering of 
photons traversing the developed emulsion. If the photographic image 
is observed over a very small area very close to the emulsion surface, 
it is not clear how the observed intensity is related to the intensity 
distribution which exposed the emulsion; in other words, the micro- 
image is not yet adequately understood. 


The approach usually followed is to model a micro-image 
scan as a continuous signal to which white gaussian noise has been added. 
Although a model postulating a photographic micro -image to be the 
result of a continuous signal having modulated a white noise "carrier" 
might be somewhat closer to the actual physics, the standard signal 
plus noise model is adopted for this study. It is emphasized that this 
approach is taken not in support of the standard model but simply 
because to establish a better model exceeds the study scope. 


For the electronic image processing system elements, 
an additive white gaussian noise model is theoretically as well as 
pragmatically acceptable. 


Subjective Image Quality 

Efforts have been made in a number of studies 53 


to define subjective assessment of image quality in a quantitative manner. 
Such techniques by definition include human variables which are not 
well defined. It is not surprising that a universal subjective measure 
has not been accepted although correlation of subject response with 
measurable parameters within the limits of specific photographic 


product uses has been shown. 
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Subjective image quality efforts are directed to 
achieving a causal relationship between measurable system element 
performance properties: frequency response (OTF), signal-to-noi.se 
ratio, and the ability of the user to make subjective judgments (usually 
in the form of detection /recognition decisions) on the output product. 
The motivation of such efforts is consistent with the objective of this 
study. However, the estimation error is the result not only of an 
estimation error introduced by the image processing system but also 
of an error introduced by that subject's own detection process. There 
is no way to separate the two unless an independent value for either can 
be obtained. 


To determine quality control procedures for the NDPF 
requires determination of the estimation error due to the system 
itself. Thus, subjective image quality measures are useless unless 
a valid model of the human subject "receiver" properties and a 
statement of its optimality were to be available. As this model is not 
established, subjective techniques necessarily lack the predictive 
validity required for this study's purposes. 

It should be recalled that some of the performance 
measures previously point out: resolution, MTFA, etc. are 

subj ective. 
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9. 


10 . 

11 . 

12 . 

13. 


14. 


15. 


16 . 


17. 


18. 

19. 


20 . 


21 . 

22 . 
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SVS-7854, "Preliminary Specification for the Data Services 
Laboratory of the NDPF" 

SVS-7756, "Performance and Design Requirements for the ERTS 
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High-Resolution Film Recorder" 

ES-0172, "Electron Beam Recorder Controller" 

ERTS-2001, "Bulk Processing Subsystem Performance Specification" 
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Controller " 

ERTS-2002, "Multispectral Scanner Video Tape Recorder /Reproducer" 
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Wmi HUM 


ruigplou cha^i sei - iRSFwy? 


FORTRAN flODULfc* 

CRRO NO 

L 

a 

3 

4 

5 
S 
7 - 
8 
9 

20 . 

1L 

12 

13 

H 

15 

16 
!7 
18 

19 

20 
2L 
22 

23 

24 

25 
2S 
27 
20 

29 

30 

31 

32 

33 
3‘l 
35 
3S 

37 

38 

39 ‘ 

-10 

41 

42 
43 - 

44 • 

45 
45 
47 
49 
49 
30 
5L 

52 

53 

54 


U.ISTJ 

**** , coin E 111 5 >Mr> 

C, -ff.ll. STEIN 3/73 ; 

C •" 

c 

C THIS PROGRAM SIMULATES I URGE PROCESSING SYSTEMS FOR THE PURPOSE OF 
G CALCULATION OF VARIANCES IN PARAMETER E5HNA1 TOM FROM NE.ASUREIOlTS- 
C MODE ON- OUT FUT_ nTQ.tiUCTS-AS-A-FuNCT I Off 'OF ' SUDSTSTCM El CHEW 

- C* PERFORMANCE. ANY 5Y5TEH WHICH CAN BE REPRESENTED CV fill ARBITRARY 

C SEQ'JEH CE OF LINEAR CLCMENTS.NONL (NEAR CHINS. AND GOOtTtVE NOISE CAN 
C. BE STHULftTED. INPUT TARGETS GENERATED HAVE THE FOWt Y.FiX.AUiJ 
C WHERE 1HE AU) ARE TARGET PARAMETERS AND X-OTSTANCE (CURRENTLY 
C T.LE.4). THE CRANER-1W0 BOUNDS ON THE VARIANCES OF THE All) 

C ESTIMATES ARE CALCULATED. 

C. 

c 

C THE PROGRAM OPERATES IK OtlE OF THREE MODES 
C ' MOBC=I IMPLIES CALCULATION OF VARIANCES 

C HG0E*2 IMPLIES CALCULATION -OF NOISE SPECTRAL DENSITY AND PRE-WHITENING 

C FILTER SPREAD FUNCTION 

C ' M0D£*3 IMPLIES TARGET 15 PLOTTED FTTCR EACH SY5TEM ELEMENT 

C UNTTSr THE UNITS OF THE CALCUALRIED VARIANCES* THE UNITS OF THE 
C RESPECTIVE AlU SQUARED 

G 

C DEFINITIONS 

C. WPUUClUlE.O YIELDS OUTPUT TARGET PUtlCHEO OH CARDS 
C H0N1F ..HUMBER OF REPEII! IONS FOR NOISE POWER SPECTRUM CALCULATION 

C HHPD-JlHtlE NOISE POWER SPECTRAL DENSITY 

C IREFC--RAHDCIH HUMBER GENERATOR SEED 

C NA*NUMBER OF A'S OF INPUT TARGET 

C OIFR-DECIMAL FRACTIONAL CHANGE tN EACH A' FOR PARTIAL DIFFERENT CAL 

G . COMFUTRT TON _ ’ . 

C IF ANY All! REPRESENTS FOSTTION. THEN UIFA.NU5T DE CHOSEN SUCH THAT; 

G DTFR4RU J/DELTX s INTEGER 

C ' Hlll=mJM6ER OF VALUES OF RUI 
C DA IT I* INCREMENT IN VALUE OF fll U INCREASE 
C fillMWTMl VALUE OF Alt) 

C CHOOSE A131.A14J EQUAL TO INTEGER MULTIPLES C p DEITY -j 

C NX -NUMBER OF P0IN1S IN TARCET GENERATED (CURRENTLY NX.LE.330! 

C 0ELTX»D13TANCE BETWEEN POINTS IN X DIRECTION 
C N8LGCN --NUMBER OF SYSTEM BLOCr.5 -ONE ST STEH ELOlM-LTNcAR ELEMENT, 

G NONLINEAR CLOICNT, AODTT IVE NOISE ICUFREHTIY (CLOCK. LE. 201 

C OPlt.Jl-0 IMPLIES SYSTEM ELEMENT TYPE J OYPASSrD IN I'TH 

C SYSTEM 8L0CK 

C J-l REFERS TO LINEAR ELEMENT 

C J*2 refers to nonlinear element 

C JOREFERS TO AUDITIVE NOISE 

C NRUN-NUK8ER OF SYSTEM BLOCK PERFORMANCE CONFIGURATIONS 
C ICOHHRUH. no IMPLIES NO CHANGE IK JFCH DATA FOR t"JM BLOCK 
t l LEV U RUN. 1 1 -0 IMPLIES NO CHANCE IN ISKil INEffiWV for rin clock 
C INOISURUN.U--0 IHTLIES NO CHANCE TN NOISE FOR t'TH BLOCK 
G NS*HUMBER OF POINTS IN S r CN 

C SFCH*SPRfOD FUWCTION OF PARICUIAR SYSTEM BLOCK 
C SET THE NUMBER OF POINTS IN EACH SFCH-’W 
C THE X INCREMENT OF SFCN MUST ALSO EQUAL 5ELTX 

C NDtlJ*HUMBER OF POINTS TN rt-B CURVE fqr BLOCK t 
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IOW/73 


INPUT USllUG 


WJT0FU3U CHART SET - CRLSPni© 


CfiRB r«3 **** 


coin E UTS 




S3 

se 

5? 

S3 

59 

60 
61 
62 

63 

64 

65 

66 
67 
GS 

69 

70 

71 

72 

73 

74 

75 

76 

77 
73 
79 


C Dim I. Jl -INf'Ul 5EMSIT7 VALUES FOR H-0 OWE OF O'-OCf t 
C OOUt It, JJ-00TPU1 PEHS11Y VRIUC8 FOR OtUll.JJ 
C 5IGHfl*STD. DEVIflUOH OF UOtSE 

DIMENSION Nl4i.flflt-tj.Ll4l.MI H I . VllGI .0141 , VI1I2O0! 

COMMON WT.NX.HOLOCii. 15. TRIM. IB.mFq.cCLlX.OPlZO.SJ. IC0m2D.2Ol. HE 
IV 120.201 .1M01S 120. 20! .Hi4 1 . HP 1 4 1 .Y120.999I ,0YDqi4-.999l .S1G14, 4i ,SF 
1CN 120, 200) . SIC-Mft 1 20) . NOPE. IERR. CIEV 1 20, 4 ) , HOME . 

loimao . 3 oj . pout 120 , boj . no 120 j , ns 
20i com 1 HUE. 

HR11EI6, 2013) 

C INHIRUZE C0H5T HN IS AND CONTROL IHfltCES 
~ REAP 15 , tool . CND-200! NODE 
RCPP 1 5. tOOt JNFUIlCrt 
RERDI5, 1001 JM0H1E 
RERD15, toil lUtlPfl 
CO 10 121. 23, 23) .MODE 

21 HR TIE 16,201 4) MODE 
CO 10 24 

22 MR UEMf» 201 5J NODE 
CO 10 24 

23 NRllE 16. 20161 CODE 

24 CONTINUE 

wmetG.?om 

READ 15, 1010) 1REFG 
CRlL CRHORGUREFCl 


SO 

81 

82 

83 

84 
SS 
CG 
87 
63 
89 
SO 
31 

92 

93 

94 

95 

96 

97 
99 
93 
100 
101 
102 

103 

104 

105 

106 
107 
103 

109 

110 


REPO IS. 1001 JNR.Wfl 
c-i.-DtFn 

RERDI5. 1001 )1M. DETL7X 
KRllElC.2000) HR. DlFfl.NX. BcLlX 
NRllEtG, 2001! IREFC 
REAP 15. LQ04 tflBLOCK.MRlIN 
WmFlfi.2012) 

DO 13 1-1. NR 

READ 15.1002) NUl.BRllJ.Rit) 

HRlTEl6,2017JNt 1 I.PRlTJ ,R[ t) 
eui-Ritj 

13 C0NT1MJE 

Ni=nm 

IFtNR.GE.2J H.1-K12! 

IFINR.C-E.3I H3-NI3I 
IFltlfl.GE.4) N4-=N 14 1 
DO 14 t-M.NBLOCK 

14 READ (5, 1004) 10F( 1, JJ , J- 1,3/ 

DO 15 I'l.KRUH 

RERD15. 1006! 1 tCOMtl. JJ . J-I.N8C0CKJ 
REAP 15, 1006 I 1 TLEVU. JJ, J=l.h8tOCKJ 
REAP 1 5 . 1006 JUNO IS 1 1 . Jl . J- 1 . N810CK7 

15 CONTINUE 

C-S0R1 14,/(P.*DEL1X) ) 2.4343 
READ LS. 1001 J N3 
CO 99 TRUH-t.NRUU 
RRITEIC, 2002! IRON 
MR WIG. 20031 

C SET 575TEM ELEMENT PERFORMANCE PflFWHETER5/FUNCT TONS 
NU5-M3-L 
DO 20 I-l.HBLOOC 
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10/0J/73 INPUT 1 1ST TUG " nUTQFLCU CHW1 SET - COLSPW@ 

CARD NO *AM COMENT5 


111 

112 

113 

114 

115 

1 16 
117 
lie 

T19 

120 

m 

122 

123 

124 

125 

126 
'127 
120 

129 

130 

131 

132 

133 

134 

135 

136 

137 
158 

139 

140 
14L 

142 

143 

144 

145 
HG 

147 

148 

149 

150 

151 
152’ 

153 

154 

155 
15G 

157 

158 

159 
130 
151 
162 

163 

164 

165 
1GG 


iFUCONunuM, u,6a.oi go to 103 

RERD 1 5 , 1003 1 ISFCtU t . Jf , J- 1, US I 

c NcanauiE sfch 

S-lSFClUt.l I*-5FCN1 T.flSJ 1/2. 

CO 50 J=2.NH5 
50 S«S*SFCHU.Jl 
S.SpDELTX 
DO 5.1 JM.HS, 

5L SFCNII. JI*5FCHU.J!/S 

mneie, 20041 1 

wueie, 20051 iSFcmc,Ji.j-i.H5J 
io3 coNTimur 

TFuieviiRim.u.ea.oi ca to ios 

RERPlS*lCQ2JHDlt I 
NDD-NDlU 

• READ 1 5 , 1003 ) l D ill 1 1 . J) , >1 . HDD 

PER0I5. 10D3HD0UT U, JJ. >t.HDDJ 
HRTTE 16,20181 i 

WRITE 16,20201 (DlNU.J>,>l,hDD> 

NRlTEi 6. 20193 1 DOUT t E . Ji , >1 . NOD) 

2018 FORlim 1 1H0. •BLOCK'. 14.21c, -N0M.1NERR LEVEL TRANSFER CURVE' J 

20 19 FORMAT 1 LH LOUT ”, 20FS. 3 1 

2020 FORMAT 1.1 H . ‘DlN‘.20FG.3l 
105 CONT INUE 

TFllnOTSURUN.il. EO.OJ GOTO 107 
RERD t5 . 1003 IS [CUR l I J , 01 R 
- SIC.'iRUJ’-SfCiVRH jX>5q,RTtDrflJ 

WRITE 16. 2007) 1.5 TOUR ill 
107 CONTINUE 
20 CONTINUE 

C ’ SET VALUES OF THE (UU 
' TFIHODE.EQ, U HRnEtG,2010) 

DO 1 U-l.Nl 
AUI-BUJ-M U-lt*OflUJ 
IFCNR.Ea.il CO TO 10 
DO 2 T2«1.H2 
R(2i*6l2J+U2-lJ*D(U2) 

'lFtNR.E0.2l GO TO 10 
DO 3 T3-1.N3 
At3UBt3)-Ml3-lM0R13J 
IFINR.EQ.31 CO TO 10 
DO 4 (4=[.N4 

■ Rl4J«B'41-U4-US-DfU4J 

10 CONTINUE - 
DO 11 IM.Hfl 
' 11 RPTlJ^niU 

IFiHOOE.NE. 1) T5=l 
lnilOOE.NE.lt GO TO G 

C COMPUTE DIFFERENTIAL TfIPUl FOR EACH Rill 
Ml*NR»l 
DO 5 !5=l.Wl 
CO TO IS, 7,6, ’9, 121. 15 
12. CONT I HUE -v 

RPt3l-fU3l 

RPl4;--D.«lPt4J 

DO TO G 


fthk 
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10 / 0 - 1/73 


INPUT 1 1ST INC 


OJTCPLOH CHHRI SET - C4LSPftrf§> 


CURB MO 

m* - contents 

107 


9 CONTI HUE- 

ioa 


Fl°t2)-tU2J 

169 


RF13I'0W1P131 

no 


00 TO 0 

171 


8 CONT lNUE 

172 


RPUJ'RllJ 

173 


RPl2]*D*APt2) 

i74 


CO TO 6 

175 


7 CONTINUE 

178 


RP ll )-»0F8P| il 

177 


8 CONI 7NUE 

173 

c 

GENERATE TARGET RUB PR0P0CR1E THROUGH SIMULATED SYSTEM 

179 


CULL CEN 

180 


lFlliaOE.tlE.2l CO 10 G4 

- 181 


DO GTU0*2. MONTE 

182 


DO 87 W*1.NX 

1S3 


G7 YlMa.HtmifU.IIHJ 

1S4 


M«0 

165 


CO CONTINUE 

166 


N’HU 

167 


15=H 

133 


64 CONTINUE 

189 


CPU SYSTEM IHPtfiCHl 

190 


IFlMOOC.EO.i) CO TO 5 

191 


I FltlC10F.ro. 3) CO 10 4 

192 . 


lFtH.NF.,M0NTE> CO TO 60 

193 


DO 66 I5*1,M0N7E 

191 


ce cnu POSPEC 

195 


CO TO 4 

196 


5 CONTINUE 

197 

c 

COMPUTE PARTIAL DERIVATIVES HRT, THE flUl 

- 193 


CALL DIF 

199 

c 

CONFUTE INTFGRRl OVER X OF FRCH PRODUCT CONOtNAT !0t( OF PORT IRL 

200 

c 

DERIVATIVES KR1. THE Alii 1AKCN THO AT R Ui'E 

201 


CPU 1NT 

202 

' 

IFltEW.ilE.OJ CO TO 4 

203 

c 

INVERT MATRIX OP INTEGRALS 

’ 204 

c 

NRI1E MQlRtX INTO VECTOR FOR INPUT ACCEPTABLE TO HtNV 1SSPI 

205 


t.'*0 

20S 


00 40 J-I.Nfl 

207 


DO 40 NUNN 

203 


K-K*l 

209 


40 VUO'SlGlli Jl 

210 


CALL HINVlViNR.BcI.L.HU 

211 

c 

EXTRACT B1AG0NRL ELEMENTS FROM VECTOR SIOREO INVERSE MATRIX 

212 


tc-1 

213 


DO 41 1=1. lift 

214 


SIC! I, T 1 «VlK) 

215 


SIGH. ll'SICU.U+M! IPO 

2 IS 


41 

217 


TMTEU3.20QSnfitU.W.!lAl OBIGINAL PAGE IS 

218 


ilRHElG, 20091 tSlCC l* 1 J * T»l,hR) OF POOR QTTATMTYj 

219 


4 CONTINUE 

220 


3 CONTINUE 

221 


2 CONTINUE 

222 


1 CONTINUE 
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10/0-1/73 INPUT LISTING 


RUTOPLOH CHART SET - CRlSPa'{§> " , 


CARR NO * 


CONTENTS 


£23 

224- 

225' 

226 

227 

228 

229 

230 
231 .' 

232 

233 

234 

235 
238 
23 ? 

238 

239 
2-10 
24L 
2 42 
243 
244 " 
245. 
248 

247 

248 

249 

250 

251 

252 

253 

254 

255 
25S 
257 
253 

259 

260 
201 
202 
2G3 

264 

2G5 

263 

267 

268 
269- 

270 

271 

272 

273 
274- 

275 

276 

277 

278 


99 CONTINUE 
IOOI F0WIIH,FB..1J 
100: FCF.l'flTU4.2F10.4) 

1003 FC B JinmOF8.3l 

1004 FORMAT 13141 
IOOG FORMAT 14012) 

1010 FORMAT I T9) 

ion format if IP-- At- - 

2000 -FORMAT tLHQ, . I2,2)i, 'DTFA*-' ■ FG.3,41. ,T4,2l, •"KIT*-." .FS.3I 

2001 FORMATUH ,."tftBPC*‘ t 19) 

2002 Fowl ILtlli 'TRU!I»*,14J 

2003 FORMAT I LUO. ’NET I SYSTEM ELEMENT PERFORMANCE OATH") 

2004 FORMAT UHO. "BLOCK’. 14.2X. ’SPCM") 

2005 FORMAT tIH . I0F8.31 " - 

2O0G FORMAT UHO, 'GLOCK' .I4-.2X, *ClEyU’.FG.3.2X, ’CIEV2«\PG. 3,2X..'CLE«. 

r,FG.3,2X.'CLEV4»’,FG.3) 

2007- FORMATUHO, "BLOCK* , T4.2J, , STCHR»".F6.3l 
2005 FORMAT tLH ,4810.31 

2009 F0RMATUH-*-,4OX,-lEl 0.31 " 

2010 FORMAT 1//3X, "All I* ,6X". "A121 \0X. "Al3l ’.GX. "HI4l ".5X. "STOlii ’.4X. "5 
!TGt2( ’.4X, "51GL3I ' .47, "SI 014 J ' I 

2011 FORMAT UHO. "CCNST0NT5 '1 

2012 FORMAT ILHO. ’Nil; '.5N. 'Dfilll ' .5X, 'AUl "I 

2013 FORMAT ILH1, "PROGRAM SIPS’) 

2014 FORMAT UHO, "MOOE". 12,27 , "CALCULATION OF VARIANCES’! 

2015 FORMAT ILHQ, "MODE", 12. 2X, "CALCULATION OF NOISE SPECTRAL DENS IFF ANO 
1 PRE-HllttcNING FILTER SPREAD PUNCH ON "I 

2016 FORMRT UHO. "MODE*. 12. 2X, "TRACE PLOT OF TARGET PF.CIPC6AT ED THROUGH S 
1 IMULAIEO SYSTEM"! 

2017 FORMAT l LH . T4-.2F10.3J 
. GO TO 201 

200 CALL- EFPLOT 
STOP 
END " 

SVBROUT THE GEN 

C ‘ THIS SUBROUTINE GENERATES ONE OF FOUR P05SIDLE' TARGETS 
C NAM IMPLIES CONSTANT LEVEL OF All 1 

C NR-2 IMPLIES TWO DELTA FCNS Oil BACKGROUND LEVEL AUl SEPEAA1ED BY 

C DISTANCE A12) 

C NA=3 IMPLIES STEP AT At3) BETIlEEH LEVELS AUl AHO AI2J 

C M=4 IMPLIES GAR OF LEVEL AI21 ON BACKGROUND OF LEVEL Rill HlTH 

C HIBTH Rl.il AND CENTER POSITION Al+I 

COMMON NA.NK.MGLOCn", 15. IRON, IB, DIFR, DEUX, OPi20,3i, IC0NI20, 20) , 1LE 
IVI20, 201. TN01S (20. 201. A<4-l, RPt4/,Y(20, 3391, OYOA 14, 9391,51014, 4;, SF 
1CN120.200T , 51 CMH 1201 . MO PE. TERR .CLEV 120. 4 1 . MONTE . 

1DU-II20,30J , BOUT l20.3Dl.NDt2Gl.NS 
GO TO- 1301 . 502. 303 , 304 1 . NO 

301 CONTINUE 

CO 306 D-'L.NX 

VU5.TXI«RPUI 
306 CONTINUE 
CO TO 300 

302 CONI IllUE 
NXX-.API2I/DELTX1-.S 
NXXstNJt-HXXJ/2 

DO 312 tX-l.HX 


t+tt 
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CWWLOH fHOm SET - CeL3P«@ 


inpqi listing 

cord no **** coi< tents *n* 


279 

312 Yl19, IXl'GPtil 


280 

YllS.WO'J-l. 


261 

IWflt-WiJC 


262 

rns.lM--!. 


233 

00 to 300 


234 

303 continue 


285 

lixx.flp 13! / DEllX-uJ. 5 


28(5 

DO 309 CX-1.HXX 


‘287 

YU$.lX)*fiPUI 


230 

305 com 1WE 


233' 

MXX*IBX*1 


290 

DO 303 IX'NXX.tlX 


291 

rit5, 1X1 -RP12J 


29 2 

308 CONTINUE 


293 

GO 10 300 


294 

304 continue 


295 

NXX--RPt3J,'BELlX-43.5 ■ 


29S 

Mm^Pt-n/ratx-o.s • 


297 

M2-HXX/2 


293 

DO 309 [X=1,NX 


293 

rilS, 1XJ =RP (U 


300 

309 CONI Ul LIE 


301 

NXXX*NXXX-N2 


302 

NXX*KXXX*i NXX 


303 

' oo 3io cx.=nxxx,nxx 


304 

X-115. 1X1--8P12) 


305 

310 com! WE 


303 

300 com (WET 


307 

RETURN 


30-3 

END 


309 

SUBRGU1 CUE SYSlEMlI’PUNCHJ 


310 

C THIS 508800111.8 PR0PQC91ES 1R9CE1 THROUGH SENUL0TE0 SYSTEM 


311 

COMSJM WT.NX.llBl.CJCK, IS, IRON, 16,0lFO.DEUX,0Piro,3MC0Kl2O,20l, 1LE 


312 

1Y120.20I , TH01S (20. 201,fU+J.flP(4i , Y120.999I . 0vDfU4-.999l .StGM, s l ,5? 


313 

lCNt20,200l,SlGM412Ql .MODE, EEHR,CLEY»20, 4J , Home, 


314 

lPtMt20.30J , BOUT ICO, 30 J.ND 120 1,(15 


315 

1FUI0DE.NE ,31 00 10 110 


316 

1S-0 


317 

COLL FHLPPPUSJ 


318 

fS't 


319 

i to com inuc 


320 

DO 101 E»I, NBLOCK 


321 

322 
■323 

IBM 

iFtopii, ij.eq.oj oo io 104 QRIGINAIj 

cru con OF POOR 

page is 

QUAXiES 

324 

IF-m0DE.NE.3l oa 10 1C4 


325 

COLL FNLPPPllSJ 


326 

104 CONI (HUE 


327 

lFtOPU.2J.EQ.Ol 00 10 106 


328 

COLL YiEY 


329 

iFincior.ne.3i ca io iog 


330 

COLL FflLFPPtlSJ 


331 

106 com inuE 


332 

lFlOPtl.3J.FQ.OI CO 10 103 


333 

COLL Ytiaisc 


334 

1FUI0DE.NS.3J CO 10 108 
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CORO NO 

335 

336 

337 

338 

339 

340 

341 
3*2 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 
353 

359 
3 GO 

361 

362 

363 

384 
365 

365 
367 

360 

369 

370 
37L 

372 

373 

374 

375 

376 

377 

378 
373 

380 

381 

382 

383 
3S-1 

385 

366 
387 
380 

389 

390 


INPUT LISTING ■ ’ AU10FLCII CHART SET - Cfll.SPOiig) 

**** ' contents 


CPU FliLPPPUS) 

108 GOTO TRUE 

ioi corn time 

IP l INPUUCH. NE.OJ .ftliD. tl5. EQ. IHM1TE17. 3000/ ml5, EX). IX-i.NXl 
3000 F0RH91 110P&.5J 
RETURN' 

EflD ‘ 

5'JBRCOT CUE' COM 
RUIHOO* H.J.H9ZUR0U3KT 
DIMENSION CFtg-391 

GCMTH- fffl. NX. (CLOCK, 15. IRON. 10, DIPfl.DCLT/. 0^120. 3J, ICONt20.2O). ILE 

]vm:0/.IW[5T20.20l.RT+/.RP , (4|,YU0.333J.07Dfll4.933/.SIO(4,4j,$P 
ICN 1 20, 200 > . S t Clffl 1 20) . MO DE, TERR. Cl EY 1 20, 4 ) , MONT E , 

1 DIN 120 . 30 ) . BOUT 1 20 . 30;*, NO 120 ) , !)3 
' N-N5/2 
DO 700 t«l.NX 
CFttj*Q 
DO 700 J*1,NS 
K-t-M-t-J 
IFtK.n.IJ K>1 

IFlK.Ol.NTtJ K-NX 

700 CF 1 1 ) *CF llJ ■‘•SFCN 1 1 B. J J*r l I5.KI 
DO 701 E-l. NX- 

701 Yii5.n-cFin^DaTx 
RETURN 

END 

SUBROUTINE YLEY 

C THIS SUBROUTINE N0BEFTE5 THE SIGNAL BY ft PHOTOGRAPHIC HURTER- 
C ' BRIFFTEID IH-D) CURVE 

COMMON Nfl.NX.f6S.0CiC, 15. IRON. lB,DlFq,DEnX.0Pt2O,3l. lC0m2O,2O). ILE 
lV'120,20(,TN0ISl2O.2O),Bl4/.0Pt4|,Y)20,999/,D'rDf||4,909l.SICi4 I 4|.SF 
1CN (20. 2003 . S 1 CMfl 120) . MODE. tERR . CLEV 120. 4 1 . MONTE. 

1 DEN 1 20, 303, BOUT 120,30) , NO (20). NS 
1»T6 
K--NOil> 

Df!flX“OlNU,K) 

DHIM'OINlt, II 
DO 900 EX=l.NX 
DY=-flLOC10(YU5. (XII 
IFipY.GE.CflflX) OV-BT ifljc 
1FIBY, LE.DHIN) BY* DM IN 
J*Z 

902 IFtOYjLE.OTHlt, Jil CO 10 SOI 
»I 

- CO TO 902 
901 CONTINUE 

os-iBY-mNiT.jjn/inimi, ji -tunci.jjij 
OY-BOUT U. J J) -.09*1 DOUT l I . J) -TJOUTl I , J J) 3 
900 YUS, IXUiCi.*+l-DU 
RETURN 
END 

SUBROUTINE YflOTSE 

COMMON NA.NX.NBLCCK, 15. TRUN, I0.DtFfl,DELTX.(Fl2O.3). ICON) 20, 201, !LE 
1V120, 201. 11)015120. 20). A 14-J.OPl 4) ,YI20. 999). DTDfl 14.999J ,5!C(4.4(,SF 
1CN120, 200) , SICMO 120) .MODE. TERR.CLEVI20, 4 ), MONTE, 
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10/04/79 mrui usiino riuiOFLOu chart set - CRLSPnriii 

CORD NO **** CONTENTS 


391 

392 

393 
334 
395 
39G 
397 
393 

399 

400 
.401 

402 

403 

404 

405 
40G 
407 
403 

409 

410 
4IL 

412 

413 

4]4 

415 

4IS 

417 

413 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 
420 

430 

431 
452 
433 
454 

435 

436 

437 
433 

439 

440 

441 

442 

443 

44<4 

445 

44G 


1B1NI20.30J, POUT 120. 30! .NOl.'Ol.NS 
T*1B 

00 001 tXM.HX 

ril 5 ,lX)-flI 5 .IXJ*U. 45 IOilAtU*GRNIOl I 

601 com inue 
RETURN 
END 

00ER0imt.E B1F 

CCH'iMOM liO.HX.HBLOCt:, 15. IRON, IB.U1FA.DELTX , 0°l20.3l, ICON 1 20, 201, TLE 
1V120, 201, 111013(20, 20). A(+l,APl4l ,V| 20,9991 , DYDAC 4-, 9931 ,SlCt4,4i,5F 
1CH 1 20, 200 1 , SI 01401201 , MODE, TERR .CLEV (20, 4 ) , HCMlE. 

10im20.30J.OlM 120 . 30), 110120 ), 1(5 
00 201 J-l.liq 

DCMHFRtfUj) 

DO 202 IX-I.HX 

DVDtlU,lXJ-tVl>l, IX! -YU. IX!) /BA 
202 COM! INUE 
201 CONI INUE 
RETURN 
END 

subroutine tm 

c 

C SUBROUTINE TNT PERFORHS TRAPEZOIDAL INTEGRATION OtER X CF EACH 
C PRODUCT CONDI HO? 1011 OP PARTIAL DERIVATIVES HOT. THE Alii 
C TflnEN TNO AT A TINE 

t 

COHTON NA.NX.NBLOCX, 15, I RUN, 18, DIFR, OELTX, OR 120, 31, ICONI 20,201, ILE 
17120,201 , !N0ISl2O,2O! ,RI4),FlPi4i ,7120. 9931 ,DT!)f>l4,939i ,51Ci4, 41 ,3F 
1CNI20. 2001, SlOMfl 1201. MODE. rERR.ClEY!20,4 1.IICHTE. 

1 OTNI20 . 30 J , BOUT 1 20 . 30 J . NO 120 1 , 1(5 
IERR--0 

DO 40! 1*1. HA 

DO 4 D 2 J--T.NA 

DSTC'DYBSll,iHDYD«U.]l-*-PYDAtT.N*I*DYDAiJ,m) , 

0570*9510/2. 

NXX’NW 
DO 403 1X*2,HXX 

BJlG-DSTMiYlV) l J . I Xf-*OY DAt 1 , IX! 

403 CONTINUE 

sicij. Jj--Dsiciocnx 
sicu,n*-sicii.j) 

tnstoic.Ji.ne.o.' co to 402 ORIGINAL PAGE IS 

fiRiTEie, 20081 iRi in, OP POOR QUALITY 

MUTEI6,4000IUJ 
IERS=1 , 

RETURN 

402 CONTINUE 
401 CONTINUE 
2003 RJRUfUtlH .4E10.3I 

4000 FORMAT l IH-*-, 40X. 'SlGf, 12. T2, '1*0, CALC 111 FIT TON TERMINATED") 

RETURN 

END 

SUBROUTINE FRPFPitSl 
C RUTHCIN- H. J.HAZUR0H5KT 

DtNEtlSTON XPt999J,YPtg99! 

COITON IlA.HX.NBlOt*. IS. TROH..tB, DlFR.OELTX,0Pl20.3J , ICQH120, 201. !LE 
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10/CM/73 INPUT UST1HC RUKJ e LOtl CUWT SFT - WlSPmtS’ 

CARD NO **** CONTENTS 


447 


IH20, 20J . TN0i5 120. 201 .01*1. OFI4I .Y120.990J , 070014.97-31 ,St01 4. 4j ,SF 

4-18 


10)120, 2001, Stem (201 .NODE, ICRR,ClEYl20, 4 1, MONTE, 

443 


IOUU20.30J.50U1 120.301, MM20I.H3 

450 


tniS.HE.OJ GO 10 10 

45L- 


5 KU7-0 

432 


COLL PLOMO..O..O! 

433 


U7.G 

454 


DO 6 - t**L,_3 

455 


T*1-2.G 

456- 


CULL MGRIDIC.S.T , 10. ,2.5,0.2,0.25.01 

457 


6 CfiLL tl0BtDt0.5,l,I0.,2,5,2.,2.5,IJ 

458 


10 KNT'KUT*! 

459 


tFim.CT.iSl CO TO 5 

4 GO 


. X0R*8.S-24JI0Dli:Nl-1.5J 

4G1 


VORt. 7, 7-2.G4M0DI IKHT-l 1 /5.31 

4G2 


CfiLL PLOIDUJR.YOO, -3) 

483 


00 15 UJ.HX 

434 


■ YP(lJ*2/rYU5,lJ 

4G5 


15 XP t ! )’2. 54-t l , 0-T J /l)5t 

438 


CfiLL L1NElYP,XP,frt,lJ 

467 


CfiLL PLOU-OT.-WR.-ai 

463 


RElURtl 

469 


EHO 

470. 


SUBRCIUJ DC POSFEC 

471 

C 

RUTHOR- II. J.NflJUROHSI.'T 

472 

C 

THIS SUBROUTINE CfiLCULfilES THE PRE-HHlTENUlG FILTER OHO 1HE HOtsc 

473 

C 

SPECIRfiL DEHSllY 

474 


DINEU5I0H Fl200!,PSU10l,PHt ITCH .flCFl 110 l.flCFfl 12201 ,P5fi( 1101 

475 


CONHOfT Nfi.HA.HBL QCK. 15, (RUN. IB. OIFA.UELTX, OP 120,31, ! COM 1 20. 201. !LC 

478 


!V 120, 201 . TDOt 5 120, 201.(1 “LC. RPiaj .Y 120,9991 , 010014,9991 ,5!Gt4,4| ,5F 

477 


tCHt20 . 200 J . S l CMO) 20 J . NO PE. t ERR . CLEV 120. 4 ) , M0H1E . 

473 


1 D INI 20. 30 ) . POUT 120,30), HD 120I.H3 

473 


1FU5.HE.1I CO 10 5 

480 


NRIJEI6. 20201 M0H1E 

481 

2020 FOPnmuMQ,-t4ome» , ,-]4i - 

432 


CALL CLEfiRlFU J,F5flU10J) 

483 


NNPtMJ.O 

484 


55-0.0 

4S5 


5 DX*DEL7X 

48G 


WUJ.O 

487 


DO 10 1*1. NX 

4S3 


fit j-cy i is. I j 

489 


10 W*RV*F(|J 

4 SO 


RV--fi\'7fW 

491 


NF*0 . ItfU 

492 


NN-2*«F-| 

4g3 


DO 30 J-t.HF 

494 


NJUIX-JM 

495 


SUN.0.0 

496 


00 20 1*1. HJ 

497 


I J* W -L 

493 


20 5UN*SUH-|Fin -RY)*-(FllJj-9V) 

495 


30 RCFUI*5UM/HJ 

500 


flCFfilNFJ-flCFUJ 

SOL 


DO 40 1=2.HF 

502 


1B-HF-W 
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10/04/73 
CORD HO 

503 

504 

505 

506 

507 
SOS 
509 

• . 510 
51 L 
sia 
513 
5H 

515 

516 

517 
SIS 
519 
570 

521 

522 

523 

524 

525 
525 

527 

528 

529 
530, 
551 

532 

533 

534 

535 

536 

537 
533 

539 

540 
54L 

542 

543 

544 

545 

546 

547 
543 

549 

550 

551 

552 

553 

554 

555 

556 

557 

558 


input Ltsnuc 


iWt* COiH E AT S 

IU=NF=I -] 

FfCMCOSll. 5703*11 -L) .'NFI 1**2 
RCFflUUl’RCFU JAFAC 
40 ncmuGi-ncFt nitric 

CRLl FOURTH l flCFfl, NN, 03 ,fi, PH, OF, Nit. 01 

FC-l./FSUl 

SUM-0. 

00 50 l=WNH 
PSUl-*FO*PSt!l 
psaui'Psai 1 1 *p$m 

50 SUM'SUtKPSl If 
SUN‘DF*tSUit-.5) 

S3=554flCFaj 

IIHP0-HU p P4(iCFUl/l2.4SUiij 
IF US. HE. MONTEl RETURN 
00 510 C-J.tiH 

510 psiu-Psom/MOine 
hnpo»impo/uowe 

SS-SS/MCHTE 

HRTIC10.2017) 

HRTTEtG. 20051 lPSllJ, 1=1. HH1 
HR VIE 16.20181 HI iPD, 55 
HR1TEI7, 101 IUlMPD 
DO SI 1=1. NH 

51 FS 1 1 1=1. /SORT IPS 11 1 J 
COLE CLEAR IPHlll.PHlUHJI 

CSC l FORDID IPS.PH.HH.DF.F.liF, 311,01 
. Ft=l./F[t|F/2t 
DO 52 1=1. NF 

52 Fni=FC-vrlU 
llRUElG, 20191 

HRITFlG, 1003) |E| l! , t=i, W) 

HR1TE17, 10031 {Fitl.C=l.UFl 
1003 FORJ!OmCF8.3l 
1011 FOAHRT IF10.31 
2005 FORWTlLri .10F6.5I 
'2017 FORmiUHO. 'hORUfil (2E0 PONER SPECTRUM' I 

2018 FORliqilUIO. •MliPD= , ,E!6.5.2!(. ‘51CM9 SQWEH=‘ .E1G.5I 

2019 rORt-miLHO. 'PRE-HHITEHING SPREAD RJIOtOH’J 
RETURN 

END 

SUBROUTINE FOURIRl DATA. HOSTS. Ox. T RU. PH C* DC LUU.NHRSH.llOPT I 
C AUTHOR- It. J.M42UR0H3KT 

OUtEIlS ION DRTfltl J . TfiO 1 1 ) , PH11 1 1 , F 19501 . 0 19501 . 0 i95GJ ,D 1950) , 
* VCQ5t57G) , V511U51GI 
1F1H1101R/2.EQ. INDRTRH1/2J 00 TO 8 
M(Wm=HD9TR*! 

DRIAUIDAIA) -DATA IUDAW- 1 1 
0 NHR5M»NDfnfl/2 
Fill = DSTfMMnRRHl 
Gill » 0.0 

FIHHARM-U « CRT At MOAT ft I 
GlNHRRM*-! J = 0.0 
DO 1 I=2.HllflRll 
MKP = NHARIUI -1 
NHM = HHflRII-1 =J 




3AL PAGE IS 
OB QUALITY 
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RUTOFLOU CHUR! SCT - CRL3Pnt<§ 


JO/C-i/ 73 
CURB NO 

539 

550 

551 

562 

563 

564 
565'' 

566 

567 
366- 

. 563 
570 
57 L 

572 

573 

574 

575 

576 

577 
573 
573 
5S0 
581 

532 

533 

584 

535 

566 

537 

553 

589 

590 

591 

592 

593 

534 

595 

596 

597 

598 

599 
COO 
601 
602 

603 

604 

605 
60S 

607 

608 

609 

610 
GIL 
612 

613 

614 


INPUT L IST1NG 

**** CONTENTS - 


Ftl! - DRTRtMHPI-* PRlRUlHHl 

1 GUI * DRTflllJHP/- PRIRlNHIU 
HI 10 "■ HHHRI1 

SUMO -- 0.0 
NIIFI -i NURRIM 
00 2 M.NtPl 

2 SUMO » SUMO *■ Fil) 

- RUI- SUMO/HWJ' 

NHCS V 0IHRRRU3I/2 
DO 3 K»l, HHCS 

VCOS (111 -COS 13, ]4 159VUC-1/ /UNO/ 

3- VSTH IK i-57tt 13. 14 159/MK- 1 J /UNO I 
BlteiFU - 0.0 
GUI ’ 0.0 
DO 5 I-2.NHFMM 
SUN! -0.0 
SUM2 ’0.0 
DO 4. J’l.MHPt 
HD Ml-lJJ-lJ-ll/iJHtHaRMl 
L ’ U-lmJ-li-O HtMtMRRM 
R5TGN - 2/4. /I-MRHM +1 
GO 10 t31.32.33.34f.R5 tCN 

31 SUM! ’ SllMt*RJHVCO$H.*ll 
SUM2 -» SUM2*GU1*VS1HIL*1I 
GO 10 4 

32 L «■ NHfiRB-l 

SOM l - SOMl *F! jm-VCOj iL*t ) 

SUM2 - SUM2*BIJJ*VStl11l«i) 

GO 10 4 ‘ 

33 L = L-UH3RH 

SUIti SUMl-FlJj+VCOStlMJ 
SUM2 ’ SUM2-CU)*VSt«tl*ll 
GO 10 4 

34 L ’ 2*NHRR|i-L 

SOMl ’ 5UMl*FlJJ*¥C0Sll-ll 
SUN2 ’ SUN2-0lJ/*VSIN(L*l/ 

4 CONTINUE 

mil ’ SUM AM 
GUI •» SUM2/HU0 

5 CONTINUE 
sum ’O.o 

DO 6 /.’I. NHF) 

6 SUN3 * SUW^FlKI-vi-L. / 

RlNHPl/ ’ 5UM3/12./HU0I 

del mi ’ i. / in. amt w»nn- i*dxj 

TFtNOPT.EQ. 1/ CO 10 60 
DO 7 M=UntPl 

TRltlHl - 0,5«QsmfltliH42-BlMl4+2J 

(Ftemi.EQ.b.) co to st 

GO PHIlMI - R1/>N2tOII<I.Rltt|J*3. 1415927*11. -SIGN! 1. ,BIM1 11 
GO 10 7 

01 Pimm - U . -S t GNU . , H IMI 1 1 + 3 . M 15927 /2. 

7 CONI t HUE 
RETURN 

80 DO 85 H-l.NHPt 
TflUlMI-fUMI 
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(WTOFLOH CHONt SET - Cf)L5Pfi!(3' 


- 10/0-1/73 

cfra no 

615 

CIS 

617 

618 

619 

620 
621 
622 

623 

624 

625 

626 
627 
623 

629 

630 

631 

632 

633 

634 

635 

636 
G37 
636 
639 
6*0 

641 

642 

643 

644 

645 
645 
647 
643 

649 

650 

651 

652 

653 

654 

655 
65G 
G57 
653 

659 

660 
661 
663 

663 

664 
G 55 
666 
667 
669 


tnpjJi lisiihc 

•ww co«i era 5 **** 

85 WUHI ’0(MI 
RETURN 

fHO 

366600 HI. 9 FOUIMVtT W. PH I . UilRRlI. OCf. NO. FUNC . t(F UH . OCLF , KU P i 
RU1H0R- If. J.t«tZUR0U5Kl . 

Oil 16113 ION TflUlSOOl .PHU 5001. RJtCllQOO). 1(5001. GlSOOl .VC0SI2551. 
fVSlH(255l,3tC E (5011 
t-MPl ' NriRRII -*1 
HUO - WWW 
IFUTLEP.EQ.OJ GO TO 4i 
00 1 t’2,NHFl 
Ft ' t-1 • 

1 StCF(lJ=3IHi3.14i593fFt/HtiOl/(3.Wl593*Fi/HHOJ. 
stc-FdJ'i.o 
DO 4 Nl.lltIPl 

PUT! - 2.pfTfi6l(lHC03(flBSIPHt (IJ) )*StGFltl 

4 eui * 2.4rTnuai'VsunPtiiitJ i^sicfu i 
CO TO 43 

41 PO 42 t-t.NHPl 

Rill « 2,*Tf(iJlllX0SinBSlPHtlUlJ 

42 8(11 - 2.*lflUU<*SlN(P(UltII 

43 9(11 - 9(11/2.0 
NFU1I -» ZWHflRH 
DELPHI. O' 12.DJ«I0»-I)ELNU) 

(TriCS^THRRM'S *1 

TJO 5 IT - 1, (TriCS 
rn *. tt-1 

vcosiui -• cost3.-Hi533 7cfu/'--ioi ORIGINAL PAGE IS. 

5 V51HHT1 -SIN13. 141593 *T|T/HHOI OF POOR QUALM 

DO 10 N - l.NFUM U 

SUM « 0.0 
UN • -limRH *11 
NHR * IRQ51UM 

do 8 riH«i.miPt 

NNHMlMRMNll-I l 
HD * NNH/ l^llriSTJ-t) 

L » 1903 INfM-2rillPHlHnR.il) 

KSIGtl 4 24 L /'MriflriM* 1 

511 * IStGtnl.NNl 

GO TO (51.52. 53.54) .TSt CM 

51 SUM » 5Utl*flltiHMVCllSlL*li*Bllal)*- SN *VStHU*-l) 

CO TO 8 

52 L ' NHRRII-l 

SliM « SUN -filNH/riVCOStl*ll*0tftHJ* SN *V5INtl*il 
00 TO 6 

53 L * t - WARM 

SUM - SUN -RINrt)4-VC05tl. i l l-B(Nrt)* SN *V5tH(L*ll 
CO 10 8 

54 L * 2j*NH 9RH -l 

SUM «. SUI1*«l?TrtJri-VC0SIL*ll -eiNHhV S!T *VSlNll'-W 

8 COITTIMUE 
10 FUNClNl * SUM 

Reiufid 

END 
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10/04/73 KtBLF OF COHTttllS Rill) REFERENCES 

own io pncE/eox name 


POIOFLOU CKflR] SET - CR15PW 
REFERENCES ISOURCE SEQUENCE 1(0. RUD PROE7BOZ) 


PACE 1 


FORTRAN NODUIE 


CHARI UTCE - 1 NT RQ DOCTORS' C0HMEM75 


CHARI UTIE - PROCEDURES 


-t 000033 J 

3.01 

201 

1000253) 

8.26 






1000254) 

3.11 

200 

COOOOC6) 

3.05 






10000711 

3.15 

21 

1000070) 

3. 14 






l 000073 J 

3.17 

22 

1000070) 

3.14 






1000075) 

3.19 

23 

10090701 

3.14 






100007SI 

3.21 

2*1 

1000072) 

3.16 

1000074 j 

3. 18 




100003 Si 

4. 10 


(C00091) 

4.15 






1000091) 

4. IS 

13 








1000094) 

4.19 


COOOG93) 

4.17 






1000095) 

4.21 


1000094) 

• 4. [g 






>000093) 

4.23 


C0C0095) 

4.21 






100009?) 

4.24 

14 








100009 7 J 

4.24 


C000097J 

4.26 






1000099) 

4.28 


C000102J 

4.34 






1000102) • 

4.34 

15 








1000106) 

4.39 


C0D0223) 

0.26' 






1009111) 

5.05 


1000140) 

G.03 






10001161 

5.10 

50 








1000116) 

5.10 


1000116) 

5.31 






1000119) 

5.14 

51 








1000119) 

5.14 


1000119) 

5.15 






1000122) 

5.20 

103 

l 00011 L) 

5.05 






1000154) 

5,35 

105 

1000123) 

5.21 






1000159) 

6.02 

107 

1 0001 35 J 

5.36 






10091401 

6.03 

20 








1000143) 

.6.06 


1000 142) 

6. 04 






IC09J-M) 

6.07 


1000222) 

8.25 






10001471 

S.iO 


1000221 i 

6. 24 






1000150) 

6.13 


1000220) 

0.23 






1000153) 

6.16- 


10002191 

8.22 






10130154) 

6.17 

10 

1000145) 

6. 03 

I00014S) 

6.J1 

>0001511 

6.14 


100015GI 

6.19 

U 








10001561 

6. 19 


(000156) 

6.20 






1 000153) 

6.23 


C QOO 157) 

6.21 






1000162) 

6.26 


(000196) 

8.01 






>0001531 

6.27 

12 

(000162) 

6.26 






1009 16 7J 

7.01 

3 

(000162) 

6.26 






1009171! 

7103 

3 

(0C0J62J 

G. 26 • 






1009175) 

7.05 

7 

(COO 162) 

6.26 • 






10001771 

7.07 

G 

(0001 53 J 

6. 23 

1OG0LG2) 

6.2G 

1090166) 

G.23 

1000170) 

1000132) 

7.11 


(000 1S3) 

7. 14 






1000133) 

7.12 

67 








>000163) 

7.J2 


COOO 163 1 

7. 13 






1000185) 

7.16 

GO 

C0G0192) 

7.22 






1000 183) 

7. 18 

64 

COOO 180) 

7.09 






1000194) 

7.24 

6S 








>009194) 

7.24 


1000134) 

7.25 






1000196) 

8.01 

5 

(000 190 J 

7.20 






1000207) 

8.07 


(000209) 

8.11 







>000174) 


7.04 
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FRC£ 2 


10/04/73 
CARP IP PRCE-'DO-X HK 


AMOTION CHARI SET - CALSPRH 
REFERENCES (SOURCE SEQUENCE HO. AND TfiOE/QO/i 


100020 Si 

S.OS 


10002091 

8. 10 

10002091 

S.09 

40 



IOQ02HJ 

8. IS 


1000215) 

8. 17 

10002 I Gl 

8.1G 

41 



1000219) 

6.22 

4 

toootsu 

7.21 

1G00220) 

S.23 

3 



l 00022 1J 

8.24 

2 



1000222J 

8.2S 

1 



10002231 

8.26 

93 



CHARI TITLE 

- NON 

-PROCEDURAL 

STATEMENTS 



CHART nice 

- SUBROUT SHE 

GEN 


IQ0025S) 

10.01 

GEN 


1000179) 

7.08-: 

1000270) 

10,02 

301 


1.000209) 

10,01 

10002721 

10.04 



1000273) 

10.05 

1000273) 

10. OS 

306 




1000275) 

10.03 

302 


(000269) 

10.01 

10002791 

10.09 

3L2 




1000279) 

10.09 



(000279) 

10. 10 

1000284) 

10.12 

303 


(0002691 

10.01 

1000237) 

to. 15 



(000263) 

10. (6 

I000C55) 

10.16 

305 




10002311 

10.19 



(000292) 

10.20 

1000292) 

10.20 

308 




1000294) 

10.21 

304 


(000269 J 

10.01 

10002991 

10.24 



(0003001 

10,25 

10003001 

10.25 

309 




1000304) 

10.23 



(000 SOS) 

10.23 

1000305) 

10.29 

310 




1000306) 

10.30 

300 


(000274) 

10.05 

CHART TITLE - NON-PROCEDURAL 

STATEMENTS 


CHART T ULC - SUER OUT TflE 

SYSTEM INPUNCH) 


1000310) 

12.01 

5VS1EH 


(000(89) 

7.'i9-: 

1000319) 

12.05 

110 


(000315) 

13.01 

1C00321) 

12.07 



1000337) 

12.24 

10003261 

12.13 

104 


(01303221 

12.03 

1000331) 

12.18 

108 


(000327) 

12. H 

1000336) 

12.23 

108 

* 

(000332) 

12. 10 

1000337) 

12.24 

101 




1 00034 OJ • 

12.28 



(000338) 

12.25 

CHART TITLE - NON -PROCEDURAL 

STATEMENTS 



CHART TITLE 

- SUBROUT (HE CON 



10003431 

14.01 

COM 

1000323) 

12.09-: 

10003511 

14.03 


(000355) 

14.12 

1000353) 

14.05 


(000355) 

(4.(1 

(000355) 

14.08 


(000354) 

(4.00 

l 000356) 
(OQ0356I 

|4.:o 

14.10 

700 

(000355) 

(4. OS 


10001951 7.20 


' t0002S3) 10. U 


(0003241 12. It 

1090323) 12. 16 

100033d) 12.2! 


1000202) 8. O') 


1000230) 10.20 
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10/0-W73 7ADI.E OF CCJrtTEtKS AND REFEF.COCES 

CORD ID PVt/BOX NAME 


fiUUIFLOH CHARI SET - CM. SPAN 
REFERENCES (SOURCE SEQUENCE HO. AMD FftGE/80) i 


pace 3 


(00035$) 14,14 701 

100035 Si 14.14 C0003S8J 14.15 

CHART TITLE - NON-FROCE DURAL STATEMENTS 


CHART 'TITLE - SUBROUTINE 

YLEV 


1000362) 

10.0! 

YLEV 

100032?) 

-1-2; 15 

10003731 

16.03 


1000361) 

.16.14 

1000375) 

ie.ro 


(000374). 

16.04 

10003761 

16.08 


(0G0375) 

16.06 

1000377) 

1G.09 

902 

1000379) 

1G. ID 

(Q00330I 

16.il 

901 

(000377) 

16.09 

1000334) 

16.13 

900 




CHART TITLE - HSN-FROCEDURRL STATEMENTS 

CHART I CUE - SLiBRQUI IMF WCHSE 

10003921 '18.01 W01SE 1000333) 12.20-X 

1000394) 18.03 (000395) 18.04 

(0003951 18.04 C01 

CHART TITLE - NOH-PROCEDURR. STATEMENTS 


CHART TITLE - SUBROUTINE DIF 


(0004031 

20,01 

dif' 

• (OOOI98J 

8.02' 

(000404) 

20.02 


(00040$) 

20. OG 

(0004061 

20.04 


1000407J 

20. OS 

(0004071 

20.05 

202 



10004031 

20.06 

201 




CHART TITLE - NON-PROCEDURAL STATEMENTS 


CHART TITLE - SUBROUTINE TNT 


(0004131 

22.01 

INI 

(000201) 

0.03- 

1000423) 

22.03 


1000433) 

22.17 

1000424) 

22.04 


(000437) 

22. 16 

(0004261 

22.06 


1000429)' 

22.07 

1000429J 

22.07 

403 



100043 7 ) 

22. 16 

402 

(000432) 

22,09 

(0004381 

22.17 

401 




CHART TITLE - NON-PROCEDURAt STATEMENTS 


CHART .TITLE - SUBROUT INE FNIPPPUSj 


1000-44 4 J 

2+.01 

FNLPPP 

1000317) 

12.03-X 1000325) 12.12-X 10003301 12.I7-X (0003351 

10004511 

24.02 

5 

1000459) 

24. 12 

(0004551 

24. OG 


1000457) 

24.09 

1000457) 

24.08 

G 



1000458) 

24,11 

10 

(000450) 

24.01 

(000404; 

24.16 


(0004CS) 

24. 18 

(0004651 

24. J 7 

15 
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pn'j £• 4 


10/04/73 106LE OP CONTENTS Rli'J REFERENCES (TUTOFLOH CHflFlJ SET - CW.SPOH 

CORO ID PBGE/BOX IU11 IE REFERENCES (SOURCE SEQUENCE NO. null PfiCE/OG/l 

CliRRI TtlLE - MOM-FROCCtiURni SUTOIEN1S 


CHRR1 ltlLE - S'lQROUl UtE P05PEC 


1000471 J 

2G.0I 

POSPffC 

1000194) 

7.24- 

t 000485 J 

26.00 

5 

1000479) 

26.01 

1000-1661 

26.08 


. (000403) 

26, 10 

(0001691 

26.03 

10 



1 0004?4 1 

26.13 


10004*39) 

26. 1*3 

(000497) 

26.13 


10004*33) 

26. 17 

I0i0049 31 

26. iG 

-20 



l 000499 J 

26.18 

30 



(000902) 

2)3,22 


1000506) 

26.24 

tO0O5O6) 

23.23 

40 



1000911) 

26.28 


1000513) 

26. 30 

(000513) 

26.29 

50 



t 00051 S) 

27.01 


10005171 

26.32 

(0005191 

27.02 

5L0 



(000519) 

27.02 


1000519) 

27.03 

1000527) 

27.13 

51 



(000527/ 

27,13 


1000527) 

27. 14 

(0005321 

27.19 

52 



(000532) 

27.19 


(000532) 

27,20 


CHfWT HUE - NON-PRCCEDURRl STR1EHEKSS 


CURfil THLE - SUBROUTINE FOUNTFKOfllfl, NORTH. OX, TRU.PHl .DELNU.NHRWI.NOPl ) 


(000545) 

29.01 

FOURTH 

10005071 

23. 25-)! 


( 00055 1J 

29.03 

8 

1000548) 

29.01 


1000557) 

29.06 


1000500) 

2S.0S 


(000560) 

29.07 

1 




1000566) 

29.11 

2 




(000565) 

29.11 


(000505) 

29.12 


1000559) 

1000570) 

29.15 

29.16 

3 

1000570) 

29.17 

OP 1 ®^ 1, page is 
0F P -°0R QUALITY" 

(000574) 

29.20 


1000538) 

30.04 

1000577) 

29.22 


1000595) 

30.02 


1000681) 

2*3.24 

31 

1000580) 

29.23 


(000584) 

29.25 

32 

1000550) 

23.23 


1009588) 

29. S3 

53 

tOOO 550) 

29.23 


100059 2 J 

30.01 

3*4 

1000550) 

29.23 


10005951 

30.02 

4 

1000583) 

29.24 

(000587) 29.25 (000591) 29.26 

1000593) 

30.04 

5 




( 00060 1) 

30.07 

6 




(000601) 

30.07 


1000601) 

30.03 


(000606) 

.70.12 


(000611) 

30.21 


1000003) 

30.14 

50 




1000(313) 

30.15 

80 

1000604) 

30.10 


1000(314) 

30.16 


(000513) 

30. IS 


1000615) 

30.17 

35 




(009(310) 

30.20 

61 

1000607) 

30.13 


1000011) 

30.21 

7 

1000609) 

30.14 



CHRRT 1 tTLF - NON-PROCEKURRl STRIENFinS 
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10/04/73 TABLE OF C0H1 ENTS flUD REFERENCES 

CURD 10 PACE /BOX NAME 


AU10FI Oil CHARI SET - CRLS p RN 


PACE 5 


REFERENCES (SOURCE SEQUENCE HO. AND PAGE/BOl.) 


CHART HI IE - SUBROUTINE FORINVlTflU. FHI .NHRFK FET.HU.FUNC.MFUH, DEI F. KLIP) 


IOOQOI9J 

32.01 

FOR IN V 

{ 006529 J 

27. JG-X 


10000201 

32.04 


10066271 

32.06 


10000271 

32.05 

1 




10005301 - 

32.09 


(006631) 

32.11 


100053 11 

32. id 

4 




(0005331 

32.12 

41 

1 000624 J 

32.02 


100053 4 ) 

32,13 


1000635) 

32.15 


1 0005331 

32.14 

42 




10005361 

32,13 

43 

(000632) 

32.11 


100004 1J 

32.18 


(0006+3) 

32.20 


100064 3J 

32.19 

5 




1 00064 SJ 

32.22 


tOOOEGG) 

33. oe 


1000549) 

32.24 


(0006G5) 

33.04 


1000565) 

32.27 

SI 

1000654) 

32.26 


1000557) 

33.01 

52 

1000654) 

32.75 


1000500) 

33.02 

53 

(000654) 

32.26 


10006631 

33.03 

54 

(00065+1 

32.26 


10006051 

33.04 

0 

(000656) 

32.27 10006591 33.01 

10006021 35,02 

10000661 

33.06 

10 





CHART TtTLE - NON-PROCEOUPRL SIRtEHEMS 
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'0/0-1/73 


OJIOHOSI CH-ffl SC T - CflL&PWl 


PACE I 


table of dtacnosttcs 

LOCATlCDl tTlAGliOSTlt 


CARI ID 

PAGE/GOx 




t 000254) 

3.11 

unrefined 

- 

■fFPLOT 1 EXTEHNAL PEFEAEUCE 

10000701 

3.25 

UIIBCF ItlED 

- 

'GRNOKG' EXTERNAL REFERENCE 

I 00021 OJ 

8.12 

UNBCFUCD 

- 

'Nll-iV' EXTERNAL REFERENCE 

10000021 

24.03 

UUPCF1NCD 

- 

‘PLOT EXTERNAL REFERENCE 

t 0000561 

24.07 

UHOEFmCO 

- 

■HC.RtD' EXTERNRL REFERENCE 

(000-157) 

24.03 

UNPEFlieD 

- 

T.CRIO' EXTERNAL REFERENCE 

(0004621 

24.14 

UNDER H€0 

- 

'PLOT' EXTERNAL REFERENCE 

(000466) 

24.18 

UlltiEFTllED 

- 

'L IKE 1 EXTERNAL REFCP.EIICE 

1000467) 

24.21 

UNDEFINED 

- 

•PLOT' EXTERNAL REFERENCE 

1000482) 

26.04 

UNDEFINED 

- 

'CLEAR' EXTERNAL PCFEPENC5, 

1000528) 

£7.15 

UHDEFIIO 

- 

'CLEAR' EXTERNRL REFERENCE 




10/04/73 


AU1OFL0H CHART SET - CfiLSPAN 


PACE 01 


ChARl' TITIC - IN1R0BUC1GRY COMMENT 5 


THIS PROGRAM: SlMULATE$- IT!RCE'Pfi'aCFSSl.YC SYSTEMS FOR THE PURPOSE OF 


CALCULATION OF VARIANCES IN PARAMETER ESI [NAT ION FROM MEASUREMENTS 
MADE ON OUTPUT PRODUCTS AS A FUNCTION OF SUSSY$-|CM CLEMENT 
PERFORMANCE. RNY.SYSTEd IlHICH CRN BE REPRESENTED 0Y Nf( ARBITRARY 
SEQUENCE OF LI NEW- ELEMENTS, NONLINEAR GAINS. AND ADDITIVE NOISE CRN 
BE SIMULATED. INPUT 1 RfiGCT S GENERATED MOVE: THE FORM: r*Fi*.Rt|M 
WERE THE PIT? ARE TARGET PARAUE1ERS RUD ^-DISTANCE (CURRENTLY 
I.LE..4N. THE CRAI'ER-RAO BOUNDS ON THE VARIANCES OF THE PUl 
ESTIMATES WE CALCULATED, 


THE PROGRAM OPERATES 1H ONE OF THREE NODES 
MODEM IMPLIES CRLCULfiTiOH OF VARIANCES 

MODE*;: IMPLIES CFiLOJlAIlON OF NOISE 5PECTRRL DENSITT miD PRE-HH 
FILTER SPREAD FUNCTION 

M0DE=3 IMPLIES TARGET IS PLOTTED AFTER EACH SYSTEM ELEMENT 
UNtTSr THE UNITS OF THE CALCUALATEO VARIANCES* THE UNITS OF THE 
RESPECTIVE Alll SQUARED 

DEFINITIONS 

NPUNCH.IlE.O YIELDS OUTPUT TARGET FUNCHED CN CARDS 

MOUTE-flUMBER OF RfcTuT.'TIOflS FOR NOISE FQMER SPECTRUM CALCULATION • 

HNPO -WHITE NOISE POWER SFECTRRl DENSITY 

TREFG-RQNDOH HUNGER GENERATOR SEED 

NA--NUN5ER OF A‘S OF INPUT TARGET 

DIFA--BECIMAL FRACTIONAL CHANCE IN EACH A FOR PARTIAL DIFFERENT |AL 

COMPUTATION 

IF ANY mil REPRESENTS POSITION. THEN OIPA MUST EE CHOI EH SUCH 7HA 
DIFAmi I J/DELT* ' INTEGER 
Ml ti 'NUMBER OF VALUES OF A 1 1 1 
DA tl /'INCREMENT IN VALUE OF AU I INCREASE 
All) -INITIAL VALUE OF Alll 

CHOOSE A13I. Pi'll EQUAL TO INTEGER MULTIPLES OF BELT* 

NX'NUMEER OF POINTS IN TARCE1 GENERATED (CURRENTLY fN.LE.200l 
OELTMI [STANCE BETHEEN POINTS IM.jc DIRECTION 

NBLOCH*UUMOER OF SYS1CH CLOCKS — ONE SYSTEM BLOC/'UICNR ELEMENT, 
NONLINEAR ELEMENT, ADDITIVE NOISE (CURRENTLY M0I.OCK.IE. 201 
OPU.JI-O IMPLIES SYS tOT ELEMENT TYPE J BYPRSSED IN CTH 
ST3TEM BLOCK 

J'l REFERS TO L (NEAR ELENENT 
J-Z REFERS TO NONLINEAR ELENENT 
J*3 REFERS TO ADDITIVE WISE 
HRUN*NUMBER OF SYSTEM BLOCK PERFORMANCE CGNF l CURAT IONS 
ICONllRUN, 1 1*0 IMPLIES NCI CHANCE IN SFCN DATA FOR 1 1 TH BLOCK 
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10/04/73 


RUKJFLQU ClIRRT SET - CRLSPflH 


pnse o 


CHERT TITLE - INTPOO'XIOR" COMMENTS 


TLEViTRUH. I J — 0 IMPLIES HO CHANGE ,IH MOM INEfRtY FOR I'lH BLOCK 

1H013HRUH. 0-0 IIIPLICS HO CtMHCE 7H NOfSE FOR 1'Tll BLOCK 

NS*NUHBER OF POINTS IN SFCH 

SFCN'SFFERO FUHC.UOM OF PRRIOIL* SYS1EH BLOCK 

SET IMS NO! ®E R OF PC11HTS 1H EACH SFCH»HX 

THE > INCREMENT OF SFCN HOST RISO EQUAL PEL!* 

NDllJ-HUttBER OF POINTS JH H-5 CURKE FOR BLOCK 1 
mtlll.JHNfJl DENSITY VFILUL'S FOR H-D CURVE OF BLOCK [ 
COUTlT.JI'OOTPOT CEIISITY VALUES FOR 01MI.JI 
SICMH^STD. BEVIftHON OF NOISE 
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OKBteWL 
OF POOR QUALITY 



10/04/73 


mil OR. OH CHARI SET - CRLSfWi 


PACE 03 


.CHART T me - PROCEDURES 


201 


00. 26- 


I 

->l 


note ot 

-jfr * t*- •* ■A'S * * if it' 

k CONTINUE * 

***-** A- * * * * * 

I 

I 

• • - - " I 02 


/ t 

t MRtTE TO DEV /• 
/ S / . 

/ VIA FORMAT / 

/ 20i3 ' / 


I' 

1 

t 

INITIAL IFF CONSTANTS 
AND CONTROL INDUES 

I 

1 03 

/ READ FROM DEV / 

/ VIA F»fil / 

/ iom / 

/ [MO THE LIST / 

I 

1 

I NOTE 04 
. k it 4r * k A is k is 
k LIST - MODE * 

kkikkkkkkisk • 
I 

• • I 

I 

is OS 
A 

is is 

* * yes 

* END OF DATA? * 

* ‘ + 

is is 

* k 
k 

• I HO 
1 

I 

I ‘ 

[ QG - 

/•"REflrFRwToEV 7 / 

/ VIA FORMAT / 

/ 1001 / 

/ INTO THE LIST / 

* r 

’ NOTE 07 
kkk k-k kkkkkk 
*■ LIST « NFiJNCH * 

’ * k-k is i is is it k'k k 
! 

I‘ 

I 03 

^/'rehI’from'dev’/ 

/ VI R FORMAT / 

/ 1001 / 

/ iNTO THE LIST / 

r 

I HOIE 09 

AAAA-Jrj-AAAA* 

A ItST ' HOME 
kkkkkiskkiskis 

I 

I 10 

^/'rerd'froit* dev”/ 

/ VIA FORMAT / 

/ 1011 / 

/ INTO THE'LIST / 

r 

- i 


200 I It 

-v * 

1 I H 

1 I EFPLOT H 

1 1 H 

! 1 H 

i t H 

*• * 

i 

1 

I 12 

is HRLT _ * 

RET URN TO SYSTEM 


— ~ >k 

1 NOTE 13 

is kk is 4- is k k is k is 

* L 1ST - MriPD k 
kkkkkkkkkkk 
I 
1 

1 L*t 


is is 

T COMPUTED CO TO I 
1 FOR MODE t* 

Hit 

1 22 3AT I 

I 23 3.13 1 

^r"******* -- 


1 

tF OUTSIDE THE RANC-E ■ 

03. H — ’•{ 

21 i IS 


/ llRTTE TO BEY / 
/ 6 / 
2 VIA FORMAT / 

/ 20 J 4 / 

/ FROM THE 1 1ST / 

I 

- I HOTc 16 

kitkkkkkkkkk 

is LIST -- MODE k 
kkkkk'kkkkkk 
I 

I 

• • • 

I 3. 21 ! 


I I I 


24 


/ 22 / 

I 

03.14 >1 

I 17 

/TmfrialEv"'/ 

/ _ _6 

/ 'VIA FORMAT / 

/ 201 S / 

/ FROM THE LIST / 


1 

I UOIE !8' 
* LIST * MODE k 

kkkkkkkkkkk- 
I ■ 

. I 

* * * 

1 3.21 ‘. 

... 24 




/ 25 / 

I 

03.14 »| 

l 19 

awuzTo 

t 0 / 

/ VIA Ff^MAT / 

/ 2010 / 

/ FROM THE LIST / 


I 

t NOiE 20 

k is k k k k k Is is kk 

k LIST « MODE k 
kkkkkkkkkkk 
I 

0J.1G*— >r 
• 24 I -MOTE 2! 


kkki-'ti-t*i'k~t 

k CO* IT 1 1 IbE k 

kkkkkkkkkkk 

1 

I 

I 12 


/ HRHE TO DEY / 
/ 5 

/ VIA FORMAT / 

/ 2011 / 


I 

i 

I 23 


/•READ Ffifj'l DcV / 

/ s / 

/ VIA FORIE-Jl / 


/iMO I Ml! "LIST / 

I 

1 

t T43TE 2“> 

k is k k k is k fs.k k k 

* L 1ST 1 lftEFC K 

kirkkkkirkkk^s 

I 

I 

1 23 

. t -* 

II. h 

1 I GRNfF’O K 

1 1 IUEFCI H 

1 1 M 

1 t H 


... * 

• I 

1 2G 

/ READ FPCrl ?Ev” 
/ S / 

/ VIR FpPJ'Ol / 

/ ICQ1 / 


/ inia the list / 

i 

I NOTE 27 

k k k k k k k k k yis 

k LIST - HQ. DtFC - 
kkkkkkkrkki 
l 


I 

I 23 

k 

I D » I. *■ OIFA 1 

+- - 4- 

1 

I 

t 29 


/ READ FftCI DEV / 
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APPENDIX C 


SIGNAL TRACE DATA 


Notes on data in this Appendix: 

(1) Progression of target through system is down consecutive 
columns, the input target signal being shown as the upper 
left-hand corner plot. 

(2) All ordinates are calibrated from 0 to 1 arbitrary full 
scale radiance unit, R 

m 

(3) All abscissas are calibrated in units of length, the size/ 
division being given before each set of data. 

(4) The plots show the effect of each ERTS image processing 
system element (excluding noise) in the MSS nominal 
configuration defined in Section 7. 1. 

The plots on pages '214 through 233. are for three parameter edge 
targets. Record length is 100. micrometers and abscissa calibration is 
10, micrometer/box. 

The plots on pages 234 through 278 are for the four parameter 
square pulse target. Record length is 150. micrometers and abscissa 
calibration is 15. micrometers/box. Note that the record length is a bit 
short to see the entire blur of the widest (105. mm) target. 

The plots on pages 279 through 283 are for the two parameter 
pseudo delta function or "resolution" target. Record length is 15 0. micro- 
meters and abscissa calibration is 15. micrometers /box. 
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The Fortran computer program EGA determines the 
line spread and transfer function of an optical system using a digital 
representation (in density, transmission or exposure) of an edge. An . 
edge is taken to be two adjacent regions each of constant density. 


The program works in exposure space; thus, if the input 
edge array is in transmission or density space the program first con- 
verts the array to exposure space, and the H-D curve must also be 
inputted. The next major portion of the program is the input array 
smoothing. This step is optional and the effective frequency cutoff of 
the smoothing is controlled by an input constant. The smoothing is 
performed by convolving a triangle with the input record. (This is 
equivalent to multiplying the input power spectrum by a sine squared . 
function. 

The next step is the computation of the line spread function; 
first according to: 



where. AE is the exposure difference between the two sides of the edge. 


Then the optical transfer function, its modulus and phase 
are calculated as: 


where 



r 1 1 1 

= ( c % + 5I n ) 
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- 1 S ( Zk .) sin 

zrrr) 
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Finally, the line spread function is normalized by the 
distance increment: 


5(Z) = 


obigmal pageb 

fig POOR QUAtB*' 
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The following table is a list, with definitions, of the 
fortran variables used in the program. In addition, a list of the fortran 
coding and an "autoflow" flow chart. of -the- progTa'm. are included. 


N 


number of points in the input edge array 


DX 


spacing (in distance) of the points in the input edge array 


L 


IND 


ID 

E 


filter control; integer; the smaller the number, - the lesser 
the smoothing; L = 1 implies no smoothing 

controls preliminary processing of input array; integer; 

IND < 0 implies transmission input to be converted 

first to density, then exposure 

= 0 implies density input to be converted to exposure 

> 0 implies exposure input 

integer label for the input edge 

input edge data array; may be exposure, density, or trans- 
mission as a function of distance 


ES 

S 

DC, EC 

FTM 

PHI 

WF 


the E array after smoothing 
spread function array 

density, exposure for the film used; inputted if density to 
exposure conversion is required 

modulation transfer function array 

phase (of the transfer function) array 

weighting factor used in smoothing operation 
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23 
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24 
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27 
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28 
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63 
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